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by 


A.  V.  Visvanathan,  Tsai -Chen  Soong,  R.  E.  Miller,  Jr. 


The  Boeing  Company,  Seattle,  Washington 


1.  ABSTRACT 

A  classical  buckling  analysis  is  developed  for  stiffened,  flat  plates 
composed  of  a  series  of  linked  flat  plate  and  beam  elements.  Plates  are 
idealized  as  multilayered  orthotropic  elements;  structural  beads  and  lips 
are  idealized  as  beams.  The  loaded  edges  of  the  stiffened  plate  are  simply 
supported  and  the  conditions  at  the  unloaded  edges  can  be  prescribed  arbi¬ 
trarily.  The  plate  and  beam  elements  are  matched  along  their  common 
Junctions  for  displacement  continuity  and  force  equilibrium  in  an  exact 
manner.  Offsets  between  elements  are  considered  in  the  analysis. 

Buckling  under  uniaxial  compressive  load  for  plates,  sections  and  stiffened 
plates  is  investigated.  Buckling  loads  are  found  as  the  lowest  of  all 
possible  general  and  local  failure  modes  and  the  mode  shape  is  used  to 
determine  whether  buckling  is  a  local  or  general  instability.  Numerical 
correlations  with  existing  analysis  and  test  data  for  plates,  sections 
and  stiffened  plates  including  boron-reinforced  structures  are  discussed. 

In  general,  correlations  are  reasonably  good. 
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2.  NOMENCLATURE 


^S.1  *Ai2  ,A22  ,A66 

b 

B11‘B12»B22*B66 
D11»D12,D22*D66 
E11»E22,G12»G23 
Gli ,G2i’G3i,G^i 
\ 


K 

i 

^xx’^yy’^zz 

J 

fk 

KRGj 

kUG*kVG,kWG,kQG 


Wku»kv 

k8,k6,..k0 

l 


hi'hi 

m 

(m22). 

MH’M12*M22 
(ni2>i 
^n22*i 
W11 ,N22 ,N12 
*11  ’*22  ’*12 


length  of  plate,  Eq.  (5*l) 

matrix  elements,  Eq.  (4.7) 

cross  sectional  area  of  beam 

width  of  plate,  Eq.  (5.1) 

elements  of  matrix,  Eq.  (4.8) 

elements  of  matrix,  Eq.  (4.9) 

orthotropic  modulus  of  elasticity,  Eq.  (4.2) 

Eqs.  (5.10a) 

distance  of  the  kth  layer  to  the  reference  plane. 

Fig.  4.1 

distance  of  kth  layer  to  the  neutral  plane,  Fig.  4.1 
moment  of  inertia 

moment  of  inertias  of  beam  about  x,  y  and  z  axes, 

Fig.  4.2 

torsional  constant,  Eq.  (4.12) 

diagonal  matrix  for  spring  constants,  Eq.  (5.57) 

spring  constants  in  the  global  coordinate  directions, 
Eq.  (5.57) 

spring  constants,  Eqs.  (5.2)  and  (5.3) 
coefficients  of  characteristic  equation,  Eq.  (5.6) 
total  number  of  layers  of  a  laminated  plate,  Eq.  (4.3) 
Eqs.  (5.8) 

half-wave  number,  Eq.  (5.1) 

Eq.  (5.20) 
moments,  Eq.  (4.6) 

Eq.  (5.21) 

Eq.  (5.22) 

stress  resultants,  Eq.  (4.5) 

in— plane  buckling  loads,  lb/in.,  Eq.  (4.l6) 

parameter,  Eq.  (5.1) 

internal  axial  load  on  beam,  Eq.  (5*4l) 
end  load  on  beam.  Fig.  4.2 
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Vq* 

Q22 

^11*^22*  e^c* 

{H1},{R2} 

R11,R12,##,R33 


[Tdl,[Tf] 


u,v,w 


V12’V21 

w, u,v 

x. y.z 

[x1],[x2],[x3±] 

[xJltlXjUXg] 

[X^JXgU^] 


1.2.3 

[  ] 

{  > 

r  j 

(  *,x 
a,  @ 

a  ,a  .0 
x*  y #  xy 

ol 

P 


Eq.  (5.19) 

lateral  loada  on  beam  in  y  and  z  directions*  Eq.  (5.4l) 
transverse  shear  on  plate,  Eq.  (4.17a)  and  Eq.  (5*12) 
orthotropicity  constants,  Eq,  (4.2) 
stiffness  constant  for  anisotropic  plate,  layer 
displacement  coefficient  matrices,  Eqs.  (5.26)  and  (5.44) 
Eq.  (5.5) 

element  (l,  l)  of  matrix  [Q^j  ]~^,.Eq*  (4.3) 
thickness  of  layer.  Fig.  4.1 
transformation  matrix,  Eq.  (5.32),  Eq.  (5*35) 
torque  on  beam,  Eq.  4.19 

displacements  of  the  neutral  plane  of  a  plate, 

Eq.  (4.4) 

Poisson’s  ratios  of  orthotropic  plate,  Eq.  (4.2) 
displacements  constants  of  beam,  Eqs.  (5.40) 
local  coordinates.  Fig.  4.1 

matrices  in  Eqs.  (5.25),  (5.29),  (5.33)  and  (5.34) 
respectively 

matrices  in  Eqs.  (5.36),  (5.37),  (5.44)  and  (5.47) 
respectively 

matrices  in  Eqs.  (5.50),  (5.53),  (5.6l)  respectively 
matrices  in  Eqs.  (5.62) 

distances  of  offset  in  y  and  z  directions.  Fig.  5.2 

distance  to  neutral  plane,  Eq.  (4.4) 

distance  of  neutral  axis  of  laminated  plate,  Eq.  (4.3) 

coordinates  in  the  principle  directions  of  orthotro¬ 
picity,  Fig.  4.1 

rectangular  or  square  matrix 

column  matrix 

diagonal  matrix 

a(  )/&x 

wave-mode  parameter,  Eqs.  (5.1)  and  (5.18) 

in-plane  stress  components,  Eq.  (4.1) 

beam  property,  Eq.  (4.1l),  which  is  the  averaged  value 
of  (stress  times  polar  moment  of  inertia) 
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t  ,Y 

e 

o 


^  to  ^ 


L 

X 


*Yxy 


r 

* 


o>,4»,n  ,p 


Superscripts 

k 

+ 


in-plane  unit  strains,  Eq.  ( U • 1 ) 
beam  rotation  about  x-axis,  Eq.  (4.19) 
displacement  coefficient,  Eq.  (5*40) 

Eqs .  (5.4l) 

in-plane  strain  components,  Eq.  (4.1) 
warping  constant,  Eq.  (4.10) 

angle  between  global  and  local  coordinates.  Fig.  5.2 
Eqs.  (5.14)  to  (5.17) 

numbering  of  lamina  layers 

quantities  related  to  side  of  plate  at  y  =  +b/2 
quantities  related  to  side  of  plate  at  y  -  — b/2 


Subscripts 

BG 

BS 

G 

(l),(2),etc. 

i 

ij 

k 

PG 

PS 

RG 

s 


u,v,v 


quantity  belongs  to  beam  element,  global  coordinates 

quantity  belongs  to  beam  element,  local  coordinates 

refer  to  global  coordinates 

element  numbers.  Figs.  5.1*  5.3*  5.4 

numbering  of  characteristic  roots 

matrix  element  subscript 

numbering  of  lamina  layers 

quantity  belongs  to  plate  element,  global  coordinate 
quantity  belongs  to  plate  element ,  local  coordinates 
restraining  forces,  ^global  coordinates,  Eq.  (5*57) 
refer  to  offset  center  "S" 

along  directions  of  x,  y,  and  z  respectively 
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3.  INTRODUCTION 


There  have  been  numerous  publications  on  general  and  local  instabilities 
of  structure  components  under  axial  compression  that  are  made  of  flat 
plates  and  beam-like  elements.  To  cite  a  few,  Ramberg  and  Levyi*  studied 
open  section  extrusions  in  which  local  instability  was  estimated  by  buck¬ 
ling  of  flanges  taken  as  plates  with  suitable  edge  conditions  and  general 
instability  analyzed  by  treating  the  extrusion  as  a  column.  Similar 
approximations  were  used  in  Goodman  and  Boyd^  and  Goodman1  s^  studies  of 
bulb-reinforced  flanges.  For  structural  sections,  such  as  Z,  T,  channels 
and  hat-type  sections,  and  isotropic  plates  stiffened  by  such  sections, 
a  usual  practice  of  analysis  is  to  treat  them  as  an  assemblage  of  flat 
plate  elements  connected  rigidly  along  the  straight  boundaries  and  each 
element  has  the  same  sinusoidal  axial  made.  The  plate  is  usually  taken 
as  infinitely  wide  and  the  constraints  on  its  sides  are  then  neglected. 

Since  a  thin  plate  is  rather  stiff  in  the  in-plane  directions,  the 
common  Junction  between  two  plate  elements  can  be  taken  approximately 
as  simply  supported  as  far  as  lateral  displacement  is  concerned.  This 
simplification  reduces  analytical  work  considerably  and  makes  possible  some 
approximate  solutions  such  as  moment-distribution  .  However,  when  this 
simple-support  assumption  is  removed,  one  needs  not  only  to  consider  the 
lateral  displacements  at  the  Junction,  but  also  the  in-plane  motions 
thereof  induced  by  buckling.  Thus,  a  rigorous  theory  should  consider 
the  stiffened  plate  as  an  assemblage  of  a  single  or  repeated  stiffened 
panels  which  consist  of  a  series  of  flat  thin  strips,  continuously  con¬ 
nected  with  each  other  at  their  edges  and  no  restraints  assumed.  In 
this  category,  one  thinks  of  flat  plates  reinforced  by  structural  sections 
and  sandwiches  with  corrugated  core.  It  is  possible  to  derive,  in  a 
rigorous  manner,  a  unified  approach  for  a  buckling  analysis  which  need 
not  distinguish  between  the  so-called  local  modes  and  the  general  modes. 

The  lowest  load  level  that  causes  any  plate  element  to  initiate  a  buckling 
deformation  before  the  others  do  (local  modes)  or  causes  several  plate 
elements  to  have  simultaneous  deformation  (general  modes)  is  the  buckling 
load  of  the  stiffened  plate  or  section.  The  unified  approach  which  seems 
to  be  the  most  exact  at  the  present  stage  of  development  can  be  repre¬ 
sented,  for  example,  by  Wittrick’s  paper^,  and  some  analyses  of  similar 
nature  but  with  various  degrees  of  exactness  and  generality  can  be  seen 
in  Refs.  6  through  12. 

The  present  analysis  brings  the  research  on  stiffened  plates  and 
sections  another  step  forward.  Here  the  basic  building  blocks  of  the 
unified  approach,  that  is,  the  flat  plate  element  and  the  beam  element 
that  is  used  to  represent  beams ,  lips  and  beads  of  flanges,  have  been 
extended  to  unidirectional,  laminated  composites  which,  of  course,  include 
isotropic  material  as  a  particular  case.  The  theory  assumes  that  the 
orthotropic  physical  properties  of  each  layer  of  the  composites  are 
given,  i.e.  ,  E;q*  ®22*  ^12 »  v12  V21  (v21  =  ^12^22 ^®11^  * 

usual  Kirchhoff-Love  assumption  regarding  plane  section  remains  plnne 


^References  are  collected  at  the  end  of  the  text 
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is  valid.  Lips  and  Leads  of  a  flange  are  regarded  as  beams  elastically- 
attached  to  the  side  of  the  flange ,  and  the  coupling  between  axial  load 
and  the  curvature  change  is  neglected.  If  the  bead  is  composite  rein¬ 
forced,  its  physical  properties  are  calculated  in  an  approximate  manner. 
The  eccentricity  of  connections  between  plate  elements  and  between  plates 
and  beads  as  well  as  conditions  of  all  unloaded  edges  of  the  plate  ele¬ 
ments  are  considered  in  the  analysis. 

Since  in  most  cases  where  fiber  reinforced  composites  are  being  used 
as  reinforcements  in  structural  sections,  the  direction  of  the  fiber 
always  runs  in  the  same  direction  as  the  axis  of  the  section,  the  stress- 
strain  equations  used  in  the  analysis  are  restricted  to  this  type  of 
orthotroplcity.  This  simplification,  which  does  not  limit  the  generality 
of  the  present  theory,  reduces  greatly  the  complexity  of  algebraic  mani¬ 
pulations  in  the  analysis  as  well  as  in  the  computer  programs.  If  one 
needs  to  accommodate  thermal  effects  and  arbitrarily  inclined  laminated 
composites  in  the  analysis,  one  needs  only  to  use  the  appropriate  stress- 
strain  relationships  in  the  beginning  of  the  derivations.  These  equations 
are  available  in  existing  literature,  for  example.  Refs.  13  and  lU. 

In  the  following  text,  the  basic  equations  for  laminated  plates  and 
beams  were  derived  and  the  geometric  and  natural  boundary  equations 
necessary  for  enforcing  continuity  and  equilibrium  along  the  junctions 
of  plate  and  beam  elements  were  developed.  Effects  of  offsets  were 
incorporated  in  the  equations •  Elastic  instability  under  uniaxial  com¬ 
pression  on  composite  reinforced  plates,  structural  sections  reinforced 
by  composites  and  beaded  stiffeners,  and  finally,  stiffened  plates  rein¬ 
forced  by  composites  were  studied.  Buckled  shapes  from  eigenvector 
solutions  are  calculated  which  can  be  used  to  ascertain  local  and  general 
type  buckling  modes.  The  accuracy  of  the  analysis  and  the  associated 
computer  programs  were  checked  first  by  correlation  with  corresponding 
isotropic  plates,  sections  and  stiffened  plates,  then  calculations  were 
extended  to  sections  and  stiffened  plates  with  composites.  The  corres¬ 
ponding  computer  programs  are  described  in  the  following  documents: 

(Refs.  I5a,b,c)  _  BUCLAP  -  Program  Documents  for  Buckling  Analysis  of 

Laminated  Composite  Plates ;  — —  BUCLAS  -  Program  Documents  for  Buckling 
Analysis  of  Composite  Reinforced  Structural  Sections;  —  BUCLASP  - 
Program  Documents  for  Buckling  Analysis  of  Composite  Reinforced  Stiffned 
Plates.  It  is  expected  that  these  computer  programs  will  be  distributed 
to  the  public  through  COSMIC,  NASA's  computer  program  distribution  point 
at  the  University  of  Georgia. 

The  main  text  of  the  report  which  contains  the  general  theory  and 
derivations  and  correlations  with  existing  analyses  and  tests  is  a  com¬ 
plete  document.  Some  peripheral  equations  and  correlations  of  secondary 
importance,  together  with  details  of  geometry  of  test  specimens,  are. 
contained  in  the  Appendix  which,  void  of  Nomenclatures  and  accompanying 
figures,  should  be  read  together  with  the  main  text  for  a  complete  under¬ 
standing. 
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4,  BASIC  EQUATIONS 


The  first  part  of  this  section  describes  the  various  material  and 
geometric  constants  used  in  the  analysis.  Next,  the  basic  equations 
for  orthotropic  laminated  flat  plates  and  laminated  beams,  which  will 
be  used  in  the  buckling  analysis  of  section  5,  are  presented. 


Plates  composed  of  orthotropic  laminas 

In  the  following  derivations,  each  lamina  is  assumed  to  be  in  a  state 
of  plane  stress.  For  an  orthotropic  lamina,  the  stress-strain  relation¬ 
ships  are  given  by^ : 


'ak  ' 

X 


L  a 


►  s 


xy 


tr 

yn 

Q12 

CM 

4 

0 

0 

4  J 


■ 


(4.1) 


where  the  superscript  k  denotes  lamina  number  and  the  elastic  constants 
are  given  by 


Qn  =  En^1  “  v2i vi2  ^ 

^22  =  E22^1  "  V21  v12  ^ 

4  *  v21E11/(1  "  V21  v12  ^  =  V12E22^1_  v21v12^ 
4  ■  °12 


(4.2) 


The  x,  y,  and  z  axes  for  the  present  orthotropic  plate  are  assumed  to  be 
identical  with  the  principle  directions  of  the  laminates  1,  2  and  3, 
respectively  as  shown  in  Figure  4.1.  A  distance  zn  locates  the  neutral 
plane  with  respect  to  an  arbitrary  reference  plane.  This  neutral  plane 
is  determined  by  calculating  the  resultant  of  the  uniaxial  forces  in  the 
laminas  for  a  constant  and  uniform  strain  across  the  thickness. 


Thus , 

*»  ■  < X  2  ‘  3r  *k(\«  ♦  V)/(,I ». 

k=l 


11 
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k=l 


r> 

ii 

&  i-i 


(4.3) 


where  is  the  first  element  of  the  matrix  [q£j]-1,  where  ] »  being 
associated  with  kth  lamina,  is  given  in  Eq.  (4.1). 
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For  an  orthotropic  laminated  flat  plate,  the  strain-displacement  rela¬ 
tionships  are 
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where  z  *  distance  from  the  point  to  neutral  plane;  u,  v,  and  w  are  the 
displacements  of  the  neutral  plane. 


Substitution  of  Eq.  (h. It)  in  Eq.  ( U .  1 )  and  integrating  over  the  thick¬ 
ness  of  the  laminas ,  the  stress  resultants  N  and  couples  M  acting  in  the 
neutral  plane  of  the  plate  can  be  expressed  in  terms  of  neutral  plane 
displacements  u,  v  and  w  as: 
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(U.6) 


where  the  A,  B,  and  D  coefficient  matrices  define  the  overall  extensions!., 
coupling  and  bending  stiffnesses,  respectively,  of  the  laminated  ortho¬ 
tropic  flat  plate.  Figure  U.l  shows  also  the  sign  conventions. 
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The  elements  of  the  A*  B,  and  D  matrices  are  given  by-^: 


$  k 
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where  h£,  shown  in  Figure 
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(4.7) 


(4.8) 


(4.9) 


the  surface  of  the  respective  lamina  (h^  =  h^  -  zn).  It  is  to  be  noted 
that  the  B  matrix  is  responsible  for  the  coupling  between  the  membrane 
stresses  and  the  bending  of  the  neutral  plane.  For  symmetrically  laminated 
orthotropic  plates  in  which  the  mid-plane  is  the  neutral  plane,  matrix  B 
vanishes  and  coupling  does  not  exist. 


If  the  orientations  of  the  fibers  of  a  lamina  are  not  parallel  or  per¬ 
pendicular  to  the  axes  of  the  plate,  the  plate  is  anisotropic.  Then  the 
A,  B,  and  D  matrices  in  Eqs.  (4.5)  and  (4.6)  will  be  fully  populated  and 
the  quantities  (Qij)  should  be  replaced  by  their  corresponding  transformed 
quantities  (designated  as  (Qij)k  in  Ref.  l4)  in  Eqs.  (4.7)  to  (4.9).  Since 
the  present  study  involves  plates,  structural  sections  and  stiffened  plates 
with  uniaxial  orthotropic  laminas  only,  these  anisotropic  equations  are  not 
given  here.  They  can  be  found,  for  example,  in  Ref.  l4.  The  specialized 
stiffness  matrices  A,  B  and  D  (Eqs.  (4.5)  and  (4.6)),  however,  can  be  used 
to  approximate  the  behavior  of  anisotropic  plates.  For  example,  when  skewed 
laminas  are  symmetric  to  the  middle  plane  of  the  plate,  skewed  orthotropic 
laminas  can  be  analyzed  in  the  present  analysis  because  these  skewing 
constants,  such  as  A-^g,  B^g,  Dqg ,  etc.,  appearing  in  Eqs.  (4.5)  and  (4.6) 
would  be  either  zero  or  very  small  compared  with  An,  Bn  and  Dn»  etc. 
and  can  be  neglected.  This  approximation  usually  requires  the  existence  of 
several  alternating  angle  plies  in  the  laminate. 


Orthotropic  laminated  beams 

Beads  or  lips  in  structural  sections,  beanntype  boron  reinforcements 
and  Joints  with  fillet,  such  as  comers  of  extruded  structural  sections, 
can  be  idealized  as  beams  and  treated  in  the  analysis.  The  origin  of  coordi 
nates  of  its  cross  section  is  chosen,  for  convenience,  at  the  geometric 
center  of  the  section.  The  basic  material  properties  involved  are  the 
individual  lamina  constants ,  E^  and  f°r  the  ^th  lamina. 

The  overall  stiffnesses  considered  for  the  laminated  beam-type  elements 
are  1^,  I^,  E^  r,  En  a  Ip  and  J.  These  are  estimated 

in  an  approximate  manner  as  follows: 
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where  F  denotes  I  ,  1  ,  I'  or  A,,  (I  and  I  are  moments  of  inertia 

yy  z*z  D  xx  yy 

about  the  neutral  axis  of  the  beam) ,  and 


~k  k 
a  I 

P 


(4.11) 


The  above  equations  are  based  on  the  assumption  that  plane  sections  of 
the  beam  cross  section  remain  plane  during  deformation.  and  Ip  for 

each  lamina  are  calculated  with  respect  to  the  axes  chosen  at  the  geo¬ 
metric  center  of  the  overall  beam  section.  The  uniform  compressive  stress 
in  kth  lamina  due  to  the  external  axial  load,  denoted  as  ak ,  is  calculated 
on  the  assumption  that  the  axial  strain  is  the  same  in  all  laminas . 

The  overall  torsional  stiffness  of  beams  whose  sections  are  made  of 
concentric  circular  layers  or  concentric  rectangular  box-type  layers  can 
be  expressed  as 

o23 ,  -  4  /  <*-12> 

For  beams  whose  cross  section  is  rectangular  and  is  composed  of  layered 
thin  plates,  Eq.  (4.12)  would  not  he  appropriate,  since  each  lamina  deforms 
with  a  different  eccentricity  towards  the  shear  center  of  the  overall 
section.  In  lacking  an  exact  torsional  stiffness  expression  for  layered 
composite  rectangular  sections,  trie  foliowing  approximate  equation  has 
been  used:  ^ 

°23  3  ‘  '<,4  4  4I(Joverall  Mrtion> 

where  is  the  cross  sectional  area  of  the  kth  layer.  Physical  proper¬ 

ties  of  these  circular  and  rectangular  cross  section  beams  can  be  found 
in  Refs.  1,  15 9  l6  and  IT* 

Equilibrium  equations  and  boundary  conditions 

The  equilibrium  equations  for  a  flat  plate  under  in-plane  compressive 
force  resultants  ^22  shear  resultant  N-^2»  written  as: 

Hii+!!i2  =  0  (4.lU) 


(4.15) 
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+  2N123x3y  =  0 


w  =  0 


+  N__3y  +  N193x  =  0 


(4.16) 


The  corresponding  boundary  conditions,  say,  at  ends  along  y— axis,  are 

SM-io  -  3w  _  3ji.  .. 

_  TIT  1  A  TtT  A  .  _  a.  A  f  li 


(4.17a) 
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ay 

=  0 

or 

M22 

=  0 

u  = 

0 

or 

N12 

=  0 

v  = 

0 

or 

W22 

=  0 

(4.17b) 

(4.17c) 

(U.lTd) 


The  equilibrium  equations  for  the  laminated,  beams  are  derived  from  ele¬ 
mentary  beam  theory  in  which  the  axial  load  due  to  externally  applied 
force  is  acting  on  the  neutral  axis  and  no  eccentricity  is  involved.  Con¬ 
sidering  a  beam  under  torsion  plus  axial  and  lateral  loads,  the  internal 
forces  (Fig.  4.2)  at  any  section  of  the  beam  can  be  written  in  terms  of 


the  applied  loads  and  displacements  u,« 

q  =  E  ,1  (d\r/dx^)  +  P  (d2w/dx2)  (4.l8) 

^z  11  yy  b 

dTx/dx  =  E1]Lr(dVdxU  )  -( G23<J-aI^(d2e/dx2)  (4.19) 

dP/dx  =  -E1 1Ab ( d2u / dx2 )  (4.20) 

q  =  E11I  (dV/dx^)  +  P  (d2v/dx2)  (4.21) 

y  11  zz  D 


where  the  rotation  0  of  the  beam  is,  in  fact,  the  same  as  the  slope  w^ 
at  the  edge  of  the  plate. 


In  the  present  analysis ,  beam  element  will  always  be  continuously  and 
elastically  connected  to  the  edge  of  a  plate  to  form  a  stiffener.  Conse¬ 
quently,  the  boundary  conditions  applicable  to  beam  elements  are: 


qz  -  022  dTx/dx  =  -M22  dP/dx  =  «12  «y  '  W22 


(4.22  -25) 


where  the  right-hand  sides  are  given  in  Eqs.  (4.17a),  (4.5)  and  (4.6). 

The  sign  conventions  are  such  that  the  beam  is  supposedly  connected  to 
the  plate  at  y  =  b  in  Fig.  4.1. 

5.  BUCKLING  ANALYSIS 

The  equations  derived  in  section  4  will  be  applied  in  this  section  to 
the  buckling  analysis  of  (a)  orthotropic,  laminated,  rectangular  plates 
and  (b)  structural  sections  and  stiffened  plates  built-up  from  laminated 
plates  and  beams.  In  the  present  linear  theory,  the  prebuckling  deforma¬ 
tions  and  possible  initial  imperfections  are  ignored.  The  buckling  load 
is  defined  as  the  smallest  load  at  which  a  part  of  the  structure  (local 
instability)  or  the  whole  of  the  structure  (general  instability)  starts 
to  develop  out-of-plane  displacements  (w-dis placement)  and  a  state  of 
unstable  equilibrium  exists  consistent  with  the  constraints  on  the  edges 
of  the  structure. 

*  Timoshenko,  S.  and  Gere,  J.  M. ,  Theory  of  Elastic  Stability,  2nd  Ed. 

1961,  McGraw-Hill,  K.  Y. ,  p.  227 
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5.1  Orthotropic  Laminated  Plate  Under  Uniaxial  Compression 

The  present  method  requires  at  least  two  opposite  edges  of  the  rectan¬ 
gular  plate  to  be  simply  supported  so  that  variables  can  be  separated  in 
the  differential  equations  of  equilibrium.  A  set  of  displacement  functions, 
automatically  satisfying  the  boundary  conditions  along  these  simply  sup¬ 
ported  edges,  are  assumed.  The  boundary  conditions  on  the  other  two  edges 
are  enforced  in  the  buckling  formulation.  Substitution  of  the  displacement 
functions  in  the  equilibrium  Eqs.  (U.lU)  to  (4.l6)  leads  to  the  character¬ 
istic  equation  of  the  differential  equations  which  in  general  is  a  polynomial 
of  degree.  Corresponding  to  each  given  load  there  are  eight  roots  from 
tnis  characteristic  equation.  The  displacements  and  the  forces  derived 
from  them  which  are  functions  of  these  roots  are  then  used  to  enforce  the 
eight  boundary  conditions  of  the  two  remaining  edges.  This  results  in  a 
set  of  eight  homogeneous  simultaneous  equations.  A  buckling  load  corres¬ 
ponding  to  each  axial  wave  number  is  obtained  from  these  equations  by 
determining  the  minimum  value  of  the  load  for  which  the  determinant  of 
the  coefficient  matrix  becomes  zero.  The  buckling  load  is  the  smallest 
of  such  loads  from  varying  the  axial  wave  numbers.  When  all  four  edges 
of  the  plate  are  simply  supported,  the  assumed  displacement  functions 
automatically  satisfy  all  the  boundary  conditions,  and  the  buckling  load 
in  this  case  is  obtained  in  a  closed  form  from  the  equilibrium  equations. 

Plates  with  two  opposite  sides  simply  supported  and  the  other  two  sides 
elastically  restrained:  _ 

As  an  example,  assume  that  a  unidirectional  compressive  load  is 
applied  along  the  direction  of  x-axis,  on  two  parallel  simply  supported 
edges  of  width  b;  the  other  two  unloaded  sides,  of  length  a,  are  elastically 
restrained.  The  four  sides  are  x  =  0  and  a,  y  =  -b/2  and  +b/2.  Let 
k  ,  k^,  k^.  be  the  translational  stiffnesses  and  kQ  the  rotational  stiffness 
of  the  elastic  restraints.  Superscripts  +  and  -  are  used  to  denote  the 
stiffnesses  and  displacements  along  the  sides  y  =  -b/2  and  +b/2  respectively. 


Displacement  functions  assumed  are: 

w  =  y  W.  sin  a  •  e^  a  =  mux/a 

.  1 


i=0L 
8 

v  *  1  V.  sin  a  *  e 

i*l  1 

8 

u  "  ][  U.  cos  a  •  e 

i=l  1 


S  =  p^y/b 


(5.1) 


The  assumed  displacements  satisfy  the  simply  supported  edge  conditions, 
v  mx1  =  N-q  =  v  =  0,  along  x  =  0  and  a.  The  conditions  at  y  =  -b/2  and 
+b/2  are: 
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where 

fin  •  -A11(m/a)2  +  A66(Pi/b)2 

R12  =  ^A12  +  AggKWaKpj/b) 

RL3  =  Bi;L(m/a)3  -  (B12  +  2B66)(m/a)(Pi/b)2 

h21  =  _R12 

R22  =  A22(Pi/b)2  “  A66(m/a)2  (5-5) 

R23  =  (B12  +  2B66)(m/a)2(P./b)  "  B22(pi/b)3 

R31  *  "R13 
**32  =  R23 

R33  =  -N11(m/ira)2  +  D^m/a)1*  -  (2D12  +  l^Dgg)  (m/a)2(Pl/b)2  +  DggCpj/b/ 


Expanding  the  determinant  equation  of  Eq.  (5.U),  one  obtains  a  polynomial: 

K8pi  +  k6p i  +  K4pi  +  K2pi  +  Ko  =  0  (5,6) 

The  above  equation  thus  yields  eight  values  of  which  are  real  or  complex 
conjugates  and  four  roots  are  the  negative  of  the  other  four# 
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From  Eq.  (5.10,  and  in  Eq.  (5.1)  can  be  expressed  in  terms  of 


(5-7) 


w. 

1 

as : 

U.  = 
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ttL0.  W. 

2i  1 

V.  = 
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R22R11 
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L2i 

R21R12  " 

R22R11 

i  ■*  1,  2 ,  •  ■  •  i 


i  =  1 ,  2 ,  . . . ,  8 


(5.8) 


The  boundary  conditions  along  the  edges  y  =  *b/2  can  now  be  enforced.  The 
displacement  functions,  with  the  eight  pj^  values  from  Eq.  (5.6)  are  substi¬ 
tuted  into  the  eight  boundary  conditions  (Eqs.  (5.2)  and  (5.3))  with  the 
aid  of  Eqs.  (4.5)  and  (4.6).  Eliminating  U*  and  V±  through  the  use  of 
Eqs.  (5.7),  one  obtains  eight  homogeneous  simultaneous  equations  in 
(i  =  1  to  8)  of  the  form  [Dg]  (W.^}  =  0.  These  equations  are  satisfied  if: 

N  *  0  (5'9) 

where  Kl  is  the  buckling  determinant  formed  by  the  coefficients  of 
This  8x8  determinant  can  be  written  as 

exp  ( wp^/2)  LG-^i  ”  ^  )  . 

exp  (irpi/2)[G2i  +  k^p^^/Ttb)]  . 

exp  (irpi/2)[G3i  -  kuL2i/ir  ]  . 

exp  (^p./2)[G1|i  -  k^L^/irJ  . 

=  0  (5.10) 

exp  (-irp^^/2)  [Gli  +  k^/ir3]  . 

exp  (-irpi/2)[G21  -  k“(pi/irb)]  . 

exp  (-rrp^^/2)  [G3i  +  ..... 

exp  (-1rpi/2)[Glti  +  k’L^/tr]  .  i  =  1,2,3,. ..,8 


i  =  1,2,3,. ..,8 
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where 


Gli  =  -B12L2i(m/a)(pi/b)+B22Lli(pi/b)2+D12(m/a)2(pi/b) 

~D2  2 ( pi /b ) 3 -2Bgg [(m/a)(pi/b)L2i+(m/a)2Lli]+4Dg6(m/a)2(piA)) 

°2i  =  -B12(Wa)L1.+B22Lli(pi/b)+D12(m/a)2-D22(pi/b)2  ( 

G3i  =  A66[L2i(Pi/b)+Lii(m/,a)]_2B66(m/a^pi/b) 

G4i  =  -A12(m/a)L2.+A22Lli(pi/b)+B12(m/a)2-B22(pi/b)2 


Since  the  p^  values  are  functions  of  the  external  load  Ny  and  the 
particular  axial  wave  number,  a  closed  form  solution  is  not  possible  and 
the  critical  value  of  tTy  has  to  be  determined  by  an  iterative  method. 

For  an  assumed  axial  mode  m  and  for  an  assumed  value  of  N]_i,  the  pi  values 
from  Eq.  (5.6)  are  substituted  into  Eq.  (5.10).  If  the  value  of  the 
determinant  is  not  zero,  the  load  is  increased  in  steps  until  Eq.  (5.10) 
is  satisfied.  This  load  then  is  the  buckling  load  for  the  assumed  axial 
mode  m.  The  calculation  is  repeated  for  various  values  of  m  and  the 
minimum  of  all  such  values  is  the  buckling  load. 

Simplifications  are  possible  for  plates  with  zero  bending-stretching 
coupling,  i.e.  By  =  0  (Eq.  (4.8))  and,  consequently,  u  and  v  are  indepen¬ 
dent  of  w  and  Eq.  (5.10)  is  reduced  to  a  4x4  determinant. 


The  criteria  of  buckling  can  be  given  in  terms  of  critical  strain  or 
critical  load  intensity.  For  composite  reinforced  stiffened  plate,  the 
critical  load  per  unit  length  is  different  for  differnt  plate  elements 
for  uniform  axial  shortening,  while  critical  strain  is  constant  throughout 
the  structure.  The  critical  strain  ex,  assuming  that  the  applied  load  is 
along  x,  is  defined  as  the  uniform  axial  strain  ex  in  Eq.  (4.1)  at  the 
instant  of  buckling  such  that  the  orthogonal  stress  Oy  in  the  kth  lamina, 
and  all  the  other  laminas,  is  zero  while  the  sum  of  the  axial  stress  aji 
through  the  thickness  of  the  plate, is  equal  to  the  applied  uniform  line 
load  I11#  Thus,  by  inverting  Eq.  (4.1)  and  putting  and  o|y  equal  to 
zero,  one  arrives  at  an  expression  of  the  critical  strain 


e 

cr 


/(  l  (tk/S*  )) 
k=l 


(5.11) 


If  one  assumes  that  the  orthogonal  stress  resultant  N22  zero,  instead 
of  lamina  stress  oK  in  each  lamina  is  zero,  one  arrives  at  the  more  familiar 
equation  for  orthotropic  plate 


ecr  =  ex  ~  ^ll^^l  "  ^A12/,A11^ 


(5.11a) 


These  two  equations  will  produce  the  same  critical  strain  criteria  since 
prehuckling  deformations  are  ignored  in  the  present  analysis. 
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Since  derivations  are  similar,  equations  for  buckling  of  composite 
plates  with  other  boundary  conditions  are  collected  in  the  Appendix  of 
the  report  for  reference,  and  will  not  be  presented  here. 

5.2  Orthotropic  Laminated  Structural  Sections  and  Stiffened  Plates  Under 
Uniaxial  Compression 

In  the  present  analysis,  structural  sections  and  stiffened  plates  of 
uniform  cross-section  are  considered  to  be  assembled  from  orthotropic 
laminated  flat-plate  and  beam  elements . 

The  intersecting  angle  between  elements  can  be  arbitrary  and  junctions 
of  elements  and  unloaded  edges  can  be  elastically  restrained,  clamped, 
simply  supported  or  free.  However,  the  loaded  edges  of  the  elements 
have  to  be  simply  supported  so  that  variables  in  the  equilibrium  equations 
can  be  separated.  For  each  element,  the  assumed  buckling  displacement 
functions  automatically  satisfy  the  simply  supported  boundary  conditions 
along  the  loaded  edges.  Consequently,  all  elements  have  the  same  axial 
mode  (wave-length).  These  simply  supported  edges  of  a  cross-section 
where  the  load  is  applied  provide  a  line  support  for  the  plate  elements 
and  a  point  support  for  the  beam  elements. 

Substitution  of  displacement  functions  such  as  Eq.  (5*1)  into  the 
equilibrium  equations  for  a  laminated  flat  plate,  Eqs.  (4,l4)  to  (4.l6), 
leads  to  a  characteristic  equation  for  each  element  such  as  Eq.  (5*6), 
which  in  general  is  a  polynomial  of  8th  degree.  Corresponding  to  an 
assumed  uniform  axial  strain  or  an  arbitrarily  variable  load  intensity 
from  element  to  element,  but  uniform  in  each  element,  there  is  a  set  of 
roots  from  this  characteristic  equation  for  each  flat  plate  element. 

For  the  beam  elements,  the  buckling  displacements  are  simple  trigonometric 
functions  in  the  axial  (longitudinal)  coordinate.  The  displacements 
assumed  are  the  translations  (u,  v,  and  w)  and  the  rotation  ($)  about 
the  longitudinal  axis.  The  enforcement  of  the  continuity  and  equilibrium 
requirements  along  common  element  junctions  results  in  a  set  of  homo¬ 
geneous,  simultaneous  equations.  A  buckling  determinant  is  formed  from 
the  coefficient  matrix  of  these  equations.  As  usual,  the  buckling  load 
of  the  structural  section,  or  the  stiffened  plate,  is  the  minimum  load 
value  among  all  ux>des  (axial  wave  numbers)  that  make  the  determinant 
vanish. 

It  is  to  be  noted  that  these  buckling  loads  and  their  eigenvectors 
are  determined  from  a  linear  exact  buckling  analysis.  Except  for  the 
axial  wave  number,  there  is  no  need  to  stipulate  a  "mode”  in  its  usual 
sense.  Eigenvector  plots  can  be  used  to  indicate  whether  the  buckling 
deformation  is  local  or  general.  In  contrast,  the  classical  buckling 
analysis  usually  assumes  one  of  the  following  individual  buckling  modes, 
in  what  is  commonly  known  flexural  (Euler)  mode,  torsional  mode,  local 
mode,  coupling  mode,  etc.  Such  simplifying  restrictions  may  result  in 
missing  the  lowest  buckling  load. 
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Element  idealization  of  structural  sections  and  stiffened  plates 

As  a  typical  example,  consider  the  arbitrary  structural  section  shown 
in  Figure  5*1*  Here  Y,  Z  are  the  global  axes  (X-axis  coincides  with  the 
axis  of  the  section  or  the  stiffened  plate)  and  y,  z  axes  are  the  local 
coordinates  of  each  component  element.  Local  coordinates  are  chosen  at 
the  geometric  center  of  beam  elements  and  neutral  planes  of  plate  elements. 
The  dash-line  contour  in  Fig.  5.1  shows  the  boundary  of  the  cross  section 
of  the  stiffened  plate;  the  solid  line  in  the  interior  is  the  neutral 
plane  of  each  segment,  from  which  offsets  between  neighboring  elements 
are  to  be  measured.  The  position  of  the  neutral  plane  can  be  calculated 
from  the  mid-plane  and  Eq.  (1+.3). 

Flat  plate  element- forces  and  displacements  in  local  coordinates 
Consider  a  typical  cross  section  of  an  orthotropic  laminated  flat 
plate  element  as  shown  in  Fig.  U.l.  The  axis  y  is  placed  at  the  middle 
of  the  plate.  The  forces  and  displacements  along  the  edges  y  =  ±b/2 
are  initially  determined  with  reference  to  the  neutral  plane.  The  forces 
involved  are : 

Q22(=  Mg2  +  2M12  ),  M22,  N12  and  N22  (5.12) 

,y  »x 

The  displacements  involved  are : 

w,  u,  v  and  v  (5.13) 

The  overall  stiffnesses  A.  ,  B.  and  D. ,  are  evaluated  as  per  Eqs.  (4.7) 
to  (4.9).  1J  1,5 


Assuming  that  the  coupling  matrix  B.  does  not  vanish,  the  buckling 


displacements  are  stipulated  as : 
8 


ij 


w  =  £  e  W.  sine*  =  £  w,  W.  sina 
i-1  1  i=l  1  1 

8  g  8 

u  =  l  irL„.ep  W.  cosa  =  J,  n.  W.  cosa 

-i  c.1l  1  •-,11 

=1  1=1 

8  6  8 
v  =  l  irL1  ,ep  W.  sina  =  £  p.  W.  sina 

i=l  -1-1  1  i=i  1  1 


and  then,  by  differentiation, 

8  6  8 

i-1  *  *  i-1 

where,  in  these  equations. 


0  *  w  *  I  ( np . /a)e^  W.  sina  *  V  W.  sina 

j  i  i  j~<  i  i 


a  -  mirx/a 


3  *  "Pjy/a 


(5.14) 

(5.15) 

(5.16) 

(5.17) 

(5.18) 
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and  symbols  n^9  <t>i#  n^,  and  P-[  are  self-evident.  Equilibrium  equations, 

Eqs •  (4.l4)  to  (4.l6),  are  satisfied  when  p-^  are  the  roots  of  the  charac¬ 
teristic  equation,  Eq.  (5*6),  and  L]i,  L2i  are  given  by  Eqs.  (5*8).  For 
the  present  axial  compression  case,  is  the  only  buckling  load  involved. 


Using  the  above  displacements,  the  forces  shown  in  Eq.  (5*12)  are 
evaluated  through  Eqs.  (4.5)  and  (4.6),  the  results  are 


^22 


i=l 


{-B12(m/a)(p1/a)L21  + 


+  D12(m/a)2(pi/a)  -  D^Cp^a)3 


-  2Bg^[(m/a)(pi/a)L2i  + 


(m/a)2Lli] 


+  4Di:'i;(m/a)2(p./a)}W.  ff^e^sina 

DO  1  1 

8 

=  l  (<l22h  Wi  sin« 

i=l 

8 

M22  =  E  {-B12^m^a^I'2i  +  B22^pi/a^LLi 

T=1 

+  D12(m/a)2  -  U22(pi/a)2}Wi„2e6sina 

8 

*  l  (m22)i  Wi  Sina 

i=l 

8 

N12  '  l  <A66[(pi/a)L2i  +  (m/a)Lli] 

i=l 

2  6 

-  2B^^(p^/a)(m/a)  }W^tt  epcosa 

8 

=  I  (ni2>i  wi  coaa 

i=l  ^  1 


(5.19) 


(5.20) 


(5.21) 
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”  l  «-A12(m/a)L21  +  A22(pi/a)Lli  +  Bl2(m/a) 

i*l 


-  B22(pi/a)2}Wi7i2e6sina 


8 

=  I  (n22^i  Wi  Sin“ 
i=l 


(5.22) 


The  quantities  (q^^i*  (*22^*  (n12^i  and  ^n22^i*  similar  to  wi»  $i% 
and  Qi  in  Eqs.  (5.36)  to  (5.39)*  are  functions  of  the  external  load 
and  the  axial  mode  m. 


Similar  equations  corresponding  tc  a  plate  with  u  and  v  not  coupled 
with  w  (Bij  =  0)  will  not  be  presented.  These  equations,  together  with 
all  such  equations  which  correspond  to  a  particular  case  of  the  general 
case  derived  here  are  collected  in  the  Appendix  for  reference. 


Transformations  to  the  reference  plane 

The  above  forces  and  displacements,  which  take  the  neutral  plane  as 
reference,  have  to  be  transferred  to  the  line  where  int  ere  lenient  matching 
is  done. 

Let  y0  and  zQ  be  defined  as  the  offsets  from  an  element  neutral  plane 
to  the  line  where  interelement  continuity  is  to  be  met.  These  offsets 
are  measured  along  the  local  axes,  with  origin  in  the  neutral  plane,  as 
shown  in  Fig.  5.2,  where  forces  and  displacements  at  point  B(y  =  +b/2) 
are  to  be  transferred  to  point  S.  The  offsets  yQ  and  zQ  shown  here  are 
taken  as  positive,  since  they  are  in  the  positive  directions  of  the  local 
axes . 


In  Fig.  5.2,  the  axes  at  point  S  are  parallel  to  the  local  coordinate 
system  y,  z  at  point  0.  The  transformation  of  the  displacements  at  B 
to  S  is  purely  a  geometrical,  rigid-body  transfer;  elasticity  between 
these  two  points  is  neglected  since  the  distance  is  small  and  conse¬ 
quently  the  effect  would  be  of  secondary  importance.  Thus ,  we  have 
displacements  at  point  S  in  terms  of  u,  v  and  w  at  B: 


w 

s 

w  +  y 

w 

s 

^o 

*y 

0 

25 

(v  ) 

-  w 

s 

,y's 

»y 

u 

55 

u  -  z 

w  + 

s 

0 

V 

— 

N 

1 

t> 

w 

s 

c 

*y 

(5.23) 


o  ,y 


The  subscript  s  indicates  the  new  displacements  at  point  S 
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Substitution  of  Eqs.  (5. 11*)  to  (5.17)  into  Eqs.  (5.23)  yields,  for 

BiJ  4  0, 

8 

v  =  l  (ok  +  yo(J>i)Wisina 
i=l 


0  =  (w  ) 

s  *y  s 


8 

£  <J>  W. sina 

i*l 


u 

s 


8 

l  [ni  -  zo(mtr/a)ui  +  y^irp^ajr^ JWjCosa 


(5.24) 


V 

s 


8 

l  (pi  "  Vi)Wisin° 


Or  written  in  matrix  form: 


w 

r  s  x 


r 


xi 
( 4x8) 

k  m 


that  i s , 

{dPS}  *  [X1]  {R1} 

the  new  notations  are  self-evident. 

In  a  similar  manner,  the  forces  along  the  edge  y 
to  a  parallel  edge  through  S,  become 


(5.25) 


(5.26) 

*  +b/2  when  transferred 


^S22^s  =  ®22  "  ZoN12,X 
(M22^s  =  M22  +  yoQ22~  ZoN22 
<B12,b  =  N12 

(N22}s  =  N22  '  yoH12,x 


(5.27) 
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and,  from  Bqs,  (5.19)  to  (5.22),  the  above  equations  become,  for  jfc  0, 


o 

{Q22)s  =  +  zo(nl2)i(m,T/a)]WiSi 


sma 


(M, 


U 

22 }s  =  l  C(m22)i  +  yo(^22}i  "  Zo(n22)i ]Wisina 


(5.28) 


8 

(N12)g  =  (n12)iW.cosa 


8 


(N22^s  *  .E  ( ( n22 ) t  +  yo(n12)i(mir/a)  ]Wisino 


which  can  be  written  in  a  matrix  form: 


'^‘Ws’ 

<VS 

<V. 


i.e. 


< 


r 


W, 


(5.29) 


{fpg}  =  [X2]  {1^}  (5.30) 

{fps>  is  self-evident. 

Equations  (5.26)  and  (5*30)  represent  the  four  displacements  and  four 
stress  resultants,  originally  acting  at  y  =  +b/2,  now  transferred  to  the 
offset  point  S  (with  positive  offsets  yG  and  zQ)  in  terms  of  the  displace¬ 
ments  at  the  neutral  plane  of  the  plate. 


Flat  plate  element -forces  and  displacements  in  global  coordinates 

Figure  5.2  also  shows  the  neutral  plane  AB  of  a  flat  plate  element  and 
its  local  coordinates  y,  z  which  makes  a  clockwise-positive  angle  with 
the  global  coordinates  X,  Y,  Z.  The  subscript  G  refers  co  the  global 
coordinate  and  the  subscript  s  indicates  the  n^cal  coordinate  at  point  S. 
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Transformed  to  the  global  axes,  the  four  displacements  of  Eqs.  (5.25) 
become : 


G> 


l  v 


G  J 


cosip 

0 

0 

-simp 


simp 

0 

0 

cosip 


Equation  (5.31)  can  be  written  simply  as 

V  “  'V  lV 


U 


Y  J 

S 


(5.31) 


(5.32) 


In  the  following  text,  superscripts  +  and  -  will  be  used  with  various 
matrix  designations  (for  example,  d+Q,  d£G,  X*,  etc.)  to  differentiate 

the  corresponding  quantities  at  sides  y  =  +b/2  and  y  -  -b/2  respectively, 
for  the  flat  plate  element.  Thus,  when  y  in  3  of  Eq.  (5*18)  is  made  to 
be  equal  to  +b/2,  we  shall  use  [X+]  in  Eq.  (5.26)  and  when  y  =  -b/2,  [Xj 

becomes  [X^].  Consequently,  substitution  of  Eq.  (5.2 6)  into  Eq.  (5.32) 

yields,  for  B-^j  #  0, 


y  =  +b/2 

(5.33) 

y  =  -b/2 

(5.34) 

displacements, 

local  coordinate 

.  to  the  global 

coordinates  by: 

cosip 

0 

0 

simp 

' <^22<«' 

(M22^G 

0 

-1 

0 

0 

< 

<“22^. 

(N12^G 

►  — 

0 

0 

1 

0 

‘Vs 

mm22)q 

0 

-sini|> 

0 

cosip 

<“22  V 

(5.35) 


Let  [Tf]  represent  the  square  transform  matrix  in  Eq.  (5.35).  Substitution 
of  Eq.  (5.30)  into  the  right-hand  side  of  Eq.  (5.35)  and  by  putting 
y  =  +b/2  and  -b/2  respectively,  one  obtains  the  following  equations  to 
represent  the  plate  forces  in  global  coordinate: 


(fpG>  =  [T^lx+KR^  =  [xj]{  R1> 
{f?G}  ^tTf3tX2]{Bl}  = 


y  =  +b/2 

(5.36) 

y  =  -b/2 

(5.3T) 
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Equations  (5.33),  (5*3*0,  (5.36),  and  (5.37)  are  displacements  and 
forces  in  global  coordinate  at  the  edges  y  =  +b/2  and  -b/2  of  a  flat 
plate  element  in  terms  of  the  eight  amplitude  constants  to  Wg  for  that 
element. 

Beam  element  -  forces  and  displacements  in  local  coordinates 

The  beam  elements  are  idealized  as  a  line  through  their  geometric  center 
in  the  longitudinal  direction.  The  forces  and  displacements  involved  in 
interelement  equilibrium  and  continuity  are  at  first  determined  along  this 
line  and  then  transformed.  The  forces  involved  are: 

q  ,  dT  /dx  ,  dP/dx  ,  and  q  (5.38) 

Z  X  jf 

The  displacements  involved  are; 

v  ,  6  ,  u  and  v  (5*39) 

The  positive  directions  of  the  displacements  and  forces  with  respect  to 
the  local  coordinate  are  shown  in  Fig.  5.2. 

The  buckling  displacements  are  assumed  as 

w  »  W  sina  0  =*  0  sina  v  -  V  sina  u  =  U  cosa  (5.*+0) 

where  a  is  equal  to  nnrx/a  as  given  in  Eq.  (5.18).  These  displacements 
satisfy  the  simply  supported  conditions  at  x  =  0  and  x  =  a  where  the 
external  compressive  force  is  applied. 

Substitution  of  Eqs.  (5.*+0)  into  Eqs.  (*+.l8)  to  (U.2l)  yields 


\  m  [E^IyyUn/a)1*  -  Pb(miT/a)2]W  sina  =  W  ^  sina 

dT^/dx  =  [Eb^r(mn/a)1+  +  G^jCmir/a)2  -  al^(mir/a)2]0  sina  =  0  sina 


dP/dx 


^EllAb^m7r^a^2^U  cosa  =  U  S 


cosa 


=  [EliIzzi(mTT/a)  -  P^mir/a)  ]V  sina  =  V  ^  sina 


(5. hi) 


Similar  to  the  flat  plate  elements,  the  transformation  of  force  and 
displacement  at  the  geometric  axis  of  the  beam  to  the  offset  axis  can  be 
done  purely  on  geometric  considerations.  Figure  5.2  can  also  be  used  to 
show  the  offset  axis  S  and  the  geometric  axis  0  of  the  beam  with  point  0 
being  moved  to  point  B. 
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The  displacements  of  Eq.  (5-39)  when  transferred  to  the  offset  axes 
become : 


w  =  v  +  y  0 

3  ^o 


u  *  u  -  z  v  -  y  v 
s  o  tx  ^o  ,x 

v  *  v  -  z  0 
s  o 


(5.42) 


where  the  subscript  s  refers  to  the  offset  axis  at  S.  After  substitution 
from  Eq.  (5.40),  the  above  equations  become: 


u 


>  = 


sina 

0 


cosa  -yo(nrrr/a)cosa  ~zo(m7r/a)cosa 


0  sina 

ur,  written  symbolically, 

{<W  *  [x5]  {r2} 


yQsina 

sina 

0 

-z  sina 
o 


(5.43) 


(5.44) 


where  [X5]  is  a  4x4  matrix.  Similarly,  the  forces  can  be  transformed  to 
the  offset  axes  as 


<«,>,  - «, + 

(dTx/dx)s  =  dTx/dx  +  qyZo  -  qzyc 

(dP/dx)  =  dP/dx 
s 


(5.45) 


Vs  ‘Sr 


+  (d2P/dx2)y, 


The  above  equations ,  on  substitution 


CO 

t* 

o* 

-z  C0(mir/a)sino 
o  3 

1  dT/dx) 

0 

X  s 

►  rr 

( dP/dx)g 

^cosa 

V* 

-y  5_(mir/a)sino 
o  3 

or,  written  symbolically,  as 


from  Eqs.  (5.4l),  become: 


0 

?1sina 

0  * 

r u 

z0Vina 

-y0^isina 

C2sina 

i 

V 

0 

0 

0 

w 

C^sina 

0 

0 

0  ‘ 

(5.46) 
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(5. ur) 


{fBS}  “  [X6]  {R2} 

where  [X^]  is  a  4x4  matrix  and  its  elements  are  functions  of  the  external 
load  on  the  beam  element.  The  displacements  and  forces  of  a  beam 
element  with  respect  to  the  offset  line  are  thus  given  by  Eqs.  (5.43) 
and  (5.46)  respectively.  The  positive  directions  of  these  are  the  same 
as  those  indicated  in  Figure  4.2. 

Beam  element  -  forces  and  displacements  in  global  coordinates 
~  Figure  5.2  can  also  be  used  to  show  the  global  coordinates  X,  Y*  and  Z 
in  relation  to  the  local  axes  y  and  2  of  the  beam  where  the  local  axis 
of  y  makes  an  angle  ^  with  the  global  Y  axis. 


Using  subscript  G  to  indicate  the  quantities  with  respect  to  the  global 
axes*  the  four  displacements  of  Eqs.  (5.43),  after  transformation*  become 


w. 


"  V 


(5.48) 


or 


{dBG}  "  [Td]  tdBS}  <5.49) 

where  the  transformation  matrix  [T^]  is  the  same  as  in  Eq.  (5.32).  Substi¬ 
tuting  for  (dgg  )  from  EQ.  (5.44)*  the  above  equation  can  be  written  as: 

{dBG}  =  [Td]  [X5]  {R2}=  IV  {R2}  (5.50) 

where  [Xy]  is  a  4x4  matrix. 

Using  the  continuity  conditions  of  Eqs.  (4.22)  to  (4.25),  the  global 
forces  acting  on  the  beam  from  the  plate  side,  whose  positive  sign  conven¬ 
tion  is  shown  in  Fig.  5.2,  can  be  expressed  by 


*^22*0  ' 

(M22)a 

*  m  IT  I  < 

{dTx/dx)s 

**12*0 

LfJ 

(dP/dx)g 

<N22)0  J 

k  (q  )  * 

7* 

*fBG*  *  *  (5.52J 


where  the  transformation  matrix  [Tf]  is  the  same  as  that  in  Eq.  (5.35). 
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(5.53) 


Substituting  for  { f^g }  from  Eq.  (5.47) ,  the  above  equation  becomes 

{fBG}  “  [Tf]  {B2}  =  tX8]  tR2} 

where  [Xg]  is  a  4x4  matrix. 

Thus,  (<igG)  and  (f“BG)  as  given  by  Eqs.  (5.50)  and  (5.53)  respectively, 
give  the  displacements  and  forces  of  the  beam  element  at  offset  S  with 
respect  to  the  global  system. 

5.3  Equations  for  Inter-element  Displacement  Continuity  and  Force 
Equilibrium 

Figure  5*3  shows  a  typical  Joint  where  three  plate  elements  Join 
together  with  a  beam  element.  The  global  coordinates  and  local  coordinates 
of  each  element  are  shown  in  the  figure  where  dash  line  illustrates  the 
outline  of  the  section.  For  the  beam  element,  B  is  its  geometric  center; 
for  the  plate  elements,  mid-plane  can  be  taken  initially  as  the  reference 
plane  for  convenience  in  measurements,  however,  position  of  neutral  plane 
(shown  in  Fig.  5.3  as  y-axis)  should  be  calculated  since  offsets  and  inter¬ 
element  matching  of  continuity  are  being  done  based  on  neutral  axes. 

Angle  ^  is  positive  when  measured  clockwise  from  Y-axis  to  y-axis  of 
the  individual  element.  For  the  purpose  of  illustration,  let  B  be  chosen 
as  the  offset  axis  on  which  inter-element  matching  will  be  done.  Conse¬ 
quently,  the  beam  element  with  center  at  B  has  no  offset. 

The  forces  and  displacements  of  the  flat  plate  elements ,  after  trans¬ 
formed  to  the  offset  point  and  converted  to  the  global  coordinates,  and 
those  of  the  beam  elements  which  are  already  in  global  coordinates,  are 
shown  in  Table  5.1. 


Table  5.1 


Forces  and  Displacements  in  Inter-Element  Matching 
at  the  Junction  in  Fig.  5*3 


Element  No. 

Force  and  Displacement 
Expression 

Equation  No. 

(l)  ana  (4) 

tfPG* 

5.37 

t(W 

5.34 

(3) 

{fPG} 

5.36 

<<G> 

5.33 

(2) 

*fBG* 

5.53 

{dBG* 

5.50 
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As  far  as  external  loads  are  concerned,  a  junction  of  elements  will 
belong  to  one  of  the  following  cases: 

(a)  free  from  external  loads 

(b)  elastically  restrained 

(c)  clamped 

(d)  simply  supported 

Detail  calculations  of  the  first  two  cases  will  be  given  as  examples. 


(a)  Free  Junction  of  elements 


Consider  the  joint  shown  in  Fig,  5*3 *  since  there  axe  no  external  con¬ 
straints,  force  equilibrium  at  B  yields 


{fPG}(l)  +  {fBG}(2)  +  {fPG}(3)  +  {fPG}(4)  “  ° 
and  inter-element  displacement  continuity  requirements  yield 

f  r,  —  -|  f  ..  i  c  1  r  -  —  1 


{dPG}(l)  =  {dBG}(2)  =  {dPG}(3)  =  {dPG}(i|) 


(5.5»0 


(5.55) 


where  the  subscripts  1,  2,  3  and  U  with  parenthesis  indicate  the  element 
number.  These  four  equations,  Eqs.  (5.5k)  and  (5-55),  after  using  the 
appropriate  equations  as  shown  in  Table  5.1,  can  be  written  in  a  matrix 
form: 


[X8](2) 

Kl(3) 

t3SJd) 

"tX73(2) 

[V(2) 

-[X3](3) 

[V(io 


[XP(3)  -[Xi](U) 


{R2)(2) 


{R1}(3) 


{V(U) 


where  the  rectangular  matrix  is  of  the  size  16x28. 


(5.56) 


(b)  Elastically  restrained  junction  of  elements 

Let  kyG».  and  kQG  be  the  four  spring  constants  of  the  restraints 

in  the  direction  X,  Y,  Z  and  0q  (Fig.  5.3)  along  the  inter-element  junction 
line.  The  appropriate  restraints  can  be  found  as 


{fRG}  =  ^RGj  tdPG}(l) 


(5.57) 


where  {fRG>  is  the  vector  of  restraint  forces  (Q22^RG»  (m22)rg*  ^N12^RG 
(^22^RGi  ^RGi  * s  “the  diagonal  matrix  whose  diagonal  elements  are  the 
spring  constants  k^G,  k^G  and  kyG;  and  {dp^^j  are  the  displacements 
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at  the  end  of  plate  No.  (l)  transformed  to  the  offset  point  B  in  global 
coordinates.  The  choice  of  the  No.  (l)  plate  is  arbitrary  since,  through 
Eqs.  (5.55),  all  elements  have  the  same  displacements  at  that  point. 


Substitution  of  Eq.  (5*3M  into  Eq.  (5-57)  yields 

(fRG>  ■  ^Gj  1*51(1)  <V<1) 


(5.58) 


which  is  the  external  force  to  be  added  to  the  equilibrium  equation  of 
the  Joint  as  given  in  Eq.  (5*5M*  The  continuity  equation  is  the  same 
as  Eqs.  (5*55)*  Consequently,  one  needs  only  to  replace  the  matrix  element 
by  [X£](l)  +  fkRQj  [X^](1)  in  Eq.  (5-56)  for  the  elastically 

restrained  Joint. 

Boundary  conditions  along  any  unloaded  flat  plate  element 


The  unloaded  flat  plate  element  can  be  free,  elastically  restrained, 
clamped  or  simply  supported.  The  constraints  and  the  displacements  will 
be  referred  to  the  neutral  plane  of  the  plate. 

By  putting  ^  =  y0  =  z0  =  0  in  appropriate  equations,  displacements  and 
forces  at  either  end  of  a  plate  element  can  be  obtained  in  global  coordinates 


(a)  Free  edge  + 

Since  forces  at  a  free  edge  should  vanish,  one  obtains  (fpG>  =  0,  or 
{f~  }  =  0.  From  Eqs.  (5.36)  and  (5*37) »  we  have 

ru 


[xu]  •  {r1)=  o 


l\]\  •  = 


(5.59) 


*  ■  y0  =  z0  =  0 

y  =  +b/2 


y  =  z 
*  o  o 


y  =  -b/2 


(b)  Elastically  restrained  edge 

Similar  to  Eq.  (5.58) ,  let  be  the  diagonal  matrix  whose  elements 

are  the  spring  constants  along  in  neutral  axis  in  global  coordinates ,  then 
the  edge  conditions  are 


([xj]  +  [k  [X+])|-{Rl}  =  0 


=  y  =  z  =  0 
Jo  o 


(5.60) 


*  fir  [^])|-{h1)=  o 


li  =  y  =  z  =  0 
r  ■  o  o 

y  =  -b/2 
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(c)  Clamped  edge 

The  boundary  conditions  for  a  clamped  edge  are  v  =  w  y  =  0  and  two 
equations  concerning  in-plane  displacements  or  forces,  i.e.,  u  or  N12 
-  0  and  v  or  N22  ®  0.  If  u  =  0  and  v  -  0  are  chosen,  for  example,  the 
clamped  edge  conditions  will  have  the  same  expression  as  Eq.  (5.59) 
except  that  Xlj  is  replaced  by  X^.  If  u  —  0  and  N22  =  0  are  chosen,  Eqs. 
(5.1*0,  (5.15),  (5.16)  and  (5-22)  should  be  used  which  will  lead  to  a 
*4x8  matrix  designated  by  a  new  notation,  say 


0 


or 


'y  -  +b/2 


(5.61) 


( d)  Simply  supported  edge 

The  classical  simple  support  conditions  are  w  =  Mg2  -  u  =  W22  =  0. 
If  these  four  conditions  are  chosen,  Eqs.  (5.lU),  (5.20),  (5.I5;  and 
(5*22)  will  lead  to  a  4x8  matrix,  designated  by 


y  *  +b/2 


or 


(5.62) 


For  convenience  of  reference,  equations  for  inter-element  matching 
and  end  conditions  are  collected  in  Table  5*2. 


Table  5.2 

Collection  of  Flat  Plate  and  Beam  Element  Equations  for 
Inter-Element  Matching  and  End  Conditions  (Bij  j*  0) 


Force  or 
Displacement 

Equation 

Equation 

No. 

Size  of 
Matrix  [X] 

Remark 

- 

tx  (R1  ) 

5.33,5.34 

4x8 

plate 

displacement 

{fPG} 

5.36,5.37 

4x8 

plate  force 

«BG>  " 

5.50 

4x4 

beam 

displacement 

<fBG>  " 

[x8]{r2> 

5.53 

4x4 

beam  force 

r— 1 

X 

tr  h- 

1 _ 1 

w-* 

II 

O 

5.59 

4x8 

free  end 

([XU]+^RG,  [X3]){R1}=° 

5.60 

4x8 

elastically 
restrained  end 

[X^]{R1}=0  or  [X*]{R1)= 

=0  5.61 

4x8 

clamped  end 

W - - 

[x*0]{Ri}=  0 

5.62 

4x8 

simply  supported 
end 
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5.4  Equations  for  Buckling  of  Structural  Sections  and  Stiffened  Plates 


The  buckling  equation  for  the  previous  structural  section  under  uniaxial 
compressive  load  will  be  obtained  by  matching  the  neighboring  elements 
together.  Assuming  that  the  left  end  of  the  section  shown  in  Fig.  5.1 
is  simply  supported  while  the  right  end  is  clamped  as  defined  in  Eq.  ( 5 • 6l ) « 

Based  on  equations  in  Table  5*2  and  Fig.  5*1*  one  obtains  the  following 
equations  which  enforce  continuities  between  elements  and  satisfy  the 
boundary  conditions  at  the  ends : 


where  the  numbers  in  parenthesis  are  the  numbering  of  elements  to  which 
the  appropriate  quantities  belong.  The  square  matrix  in  Eq.  (5.  3)  is 
the  buckling  determinant,  of  order  56x56.  A  common  factor  of  sina  and 
cosct  (a  =  mrrx/a)  should  be  taken  out  from  the  rows  of  the  determinant. 

A  second  example  considered  is  that  of  a  hat— section  stiffened  plate 
shown  in  Fig.  5.4  where  the  stiffener  dimensions  and  their  spacing  is 
repetitive.  The  stiffened  plate  is  idealized  in  Fig.  5*4  into  line  elements 
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with  the  element  numbers  shown  in  parenthesis.  The  ipter-element  continuity 
and  equilibrium  equations,  together  with  end  equations  for  the  end  elements 
(l)  and  (29) ,  can  be  written  in  a  manner  identical  to  Eq.  (5.63).  Because 
of  repetitive  nature  of  the  stiffeners,  these  equations  can  be  considered 
to  consist  of  three  basic  parts;  namely: 

(i)  a  set  of  equations,  say  designated  as  [Ta]  {Ra>,  representing  the 
left  end  side  of  the  stiffened  plate, 

(ii)  a  second  set  of  equations,  [T^]  {R^}  representing  the  repetitive 
unit.  The  number  of  repetitions  of  this  set  of  equations  is  the  same  as 
the  number  of  the  repetitive  units  in  the  structure, 

(iii)  a  third  set  of  equations,  [Tc]  {Rc>,  representing  the  right  end 
side  of  the  stiffened  plate. 

In  Fig.  5*^»  the  two  ends  and  the  inter-element  joints  are  marked  with 
a  dot.  The  left  side  block  consists  of  six  dots  and  involves  elements  from 
(l)  to  (7).  The  first  repetitive  unit,  as  well  as  all  the  other  repetitive 
units,  contains  six  dots  and  the  continuity  equations  involve  elements  from 
(7)  to  (l4).  Assuming  that  there  are  three  such  repetitive  units  in  the 
plate,  the  right  side  block  contains  two  dots  and  involves  elements  (28) 
and  (29). 

These  three  matrices  can  be  written  as  in  the  following: 


(1)  (2)  (3)  (4)  (5)  (6)  (7) 
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(7)  (8)  (9)  (10)  (11)  (12)  (13)  (Ik) 


The  buckling  determinant  can  be  obtained  from  sub-matrices  [Ta],  [T^]  and 
[Tc]t  and  the  eigenvectors  are  the  column  matrices  (Ra),  and 

As  in  Eq«  (5*63),  a  common  factor  in  each  row  for  sina  or  cosa  can  be  taken 
out  in  the  determinants!  equation  and  the  buckling  load  tod  the  corresponding 
eigenvectors  are  functions  of  the  axial  half-wave  number  m*  The  plot  of 
the  deflected  shape  by  means  of  the  eigenvectors  are  useful  for  determining 
whether  buckling  involves,  on  a  relative  magnitude,  only  local  elements 
(local  instability)  or  all  elements  (general  instability).  Computer  sub¬ 
routines  for  finding  eigenvectors  are  commonly  available  and  such  procedure 
will  not  be  described  here. 

The  structure  shown  In  Fig.  5 • k  used  for  the  second  example  is  a  typical 
stiffened  plate  which  has  a  number  of  repeated  substructures.  The  choice 
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of  how  many  elements  go  to  the  beginning  eigenvector  {Ra}  and  how  many  go 
to  the  ending  eigenvector  {Rc }  is  arbitrary.  However,  in  the  subroutine 
(named  BLKDET )  used  in  the  program  BUCLASP  (Ref.  15c),  a  working  space  is 
provided  to  first  store  [Ta]  and  decompose  it  and  then  repeatedly  bring 
in  [Tfc]  for  operation  and  finally  [Tc].  Consequently,  core  storage  can 
be  saved  by  making  [Tc]  as  small  as  possible.  That  is  the  reason  why 
[Tc]  is  smaller  than  [Ta]  in  the  second  example.  Some  techniques  in 
handling  large  matrices  and  the  convention  for  angles  relating  the  orienta¬ 
tion  of  one  plate  element  to  another  are  described  in  program  documents. 

Ref.  15a,  b,  c,  which  will  be  distributed  by  NASA. 

The  present  method  has  been  used  to  correlate  with  other  existing  analyses 
and  tests  for  composite  plates,  structural  sections  and  stiffened  plates 
in  the  next  sections.  As  has  been  mentioned  before,  the  usual  practice  of 
estimating  the  elastic  instability  of  thin-walled  structures  is  to  guess 
possible  inodes  based  on  which  simplified  analysis  can  be  made.  Junctions 
of  internal  members  are  usually  taken  as  nodal  lines  about  which  force 
and  moment  equilibrium  are  satisfied  to  some  extent  but  their  translations 
are  ignored.  The  assumption  regarding  possible  modes  requires  expertise 
in  the  analysis  and  may  miss  a  lower-load  mode  and  the  immovable  nodal 
line  assumption  requires  non-existing  external  forces  to  hold  the  joints 
which  will  make  the  structure  stiffer  than  it  is.  However,  when  the  mode 
is  correctly  assumed,  the  usual  method  may  yield  good  agreement  with  more 
exact  analysis  such  as  the  present  one,  as  can  be  seen  from  the  next  sections. 

6.  CORRELATIONS  WITH  OTHER  ANALYTICAL  RESULTS 

In  correlation  studies  with  other  analytical  results,  the  mathematical 
models  used  in  othersf  analyses  will  be  followed,  i.e.,  whether  a  stiffener 
is  treated  as  a  beam  or  a  plate  in  the  present  analysis  will  be  consistent 
with  those  in  the  references. 


A.  Buckling  of  simply  supported  web  with  isotropic  or  orthotropic  flange 
('Ref.  16) 

In  Ref.  18,  the  buckling  of  a  supported  web  section,  one  side  of  its 
cross  section  simply  supported  and  the  other  parallel  side  free,  with 
and  without  unidirectional  composite  reinforcement  on  its  flange,  is 
studied.  Two  types  of  failure  modes  were  examined  in  Ref.  18,  one  is 
a  local  mode  and  the  other  is  a  long-wave  mode. 


The  orthotropic  constants  used^-®  are 

E  =  E  in  the  direction  of  filament  -  30.25  x  10^  psi 
x 

E  =  2.03  x  106  psi 

y 


V 

xy 

Gxy 


=  0.346 

c 

=  0.521*9  X  10 


psi 


(6.1) 


unless  otherwise  specified,  these  constants  will  be  used  in  the  rest  of 
correlations  where  ortho tropic  materials  are  used. 
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Local  mode 

When  the  dimension  of  the  section  is  such  that  the  axial  wave  length 
is  of  the  order  of  magnitude  of  the  width  of  the  web,  one  observes  the 
so-called  local  mode  and  the  results  are  given  in  Fig.  6.1  where  the 
flange  is  treated  as  a  plate.  It  can  be  seen  from  the  curves  that  when 
the  section  is  reinforced  on  the  outstanding  flange  by  composites,  the 
flange  seems  stiff  enough  to  force  a  node-line  at  the  Junction  of  the 
two  legs.  As  a  consequence ,  the  approximate  analysis  in  Ref.  18  which 
implies  such  an  assumption  agrees  very  well  with  the  present  method. 

On  the  other  hand,  when  the  flange  is  not  reinforced,  the  junction  seems 
flexible  and  may  be  displaced  and  the  approximate  analysis  yields  higher 
values  than  the  present  method. 

Long-wave  mode 

For  the  long-wave  mode,  the  same  type  of  composite  reinforced  flange- 
web  is  studied  where,  in  Ref.  18  as  well  as  in  the  present  method,  the 
flange  is  treated  as  a  beam.  The  geometry  of  the  15-inch  long  column 
is  given  in  the  following, 

b  =  1.5  in.  b  =  0.6  in.  t  =  .05  in.  t  =  3  t 

W  I  W  X  w 

For  the  given  dimensions,  one  half-wave  yields  the  lowest  buckling 
load.  The  present  method  gives  a  buckling  parameter16  k^  »  1.58  while 
the  value  in  Ref.  18  is  1.56.  The  difference  is  1 %• 

B*  Buckling  of  discretely  stiffened  isotropic  plate  with  a  single  orthotro¬ 
pic  stif f ene r  ( Ref.  loT 
Shallow  orthotropic  eccentric  stiffener 

Figure  6.2  shows  the  cross  section  of  a  discretely  stiffened  plate 
with  an  orthotropic  eccentric  stiffener.  The  stiffener  is  so  shallow 
that  in  Ref.  18  it  is  treated  as  a  beam  with  Gyz  -  0.  The  boundary 
conditions  along  the  unloaded  sides  of  the  plate  are  w  =  M22  ~  ^12  = 

N22  *  0.  In  Ref.  18,  displacements  at  the  junction  between  the  stiff¬ 
ener  and  the  two  parts  of  the  plate  were  matched  for  out-of-plane 
displacement  only.  The  results  in  Fig.  6.2  showed  excellent  agreement 
in  spite  of  the  difference  in  the  rigor  of  the  two  analytical  approaches. 


Deep  orthotropic  eccentric  stiffener 

When  the  stiffener  is  deep,  it  is  taken  as  an  orthotropic  plate  in 
Ref.  18  where  forces,  moments  and  rotations  are  matched  at  the  common 
boundary  of  the  three  plate  elements  whose  lateral  displacements  at 
that  joint  are  ignored.  In  Fig.  6.3,  two  points,  marked  as  1  and  2, 
are  obtained  by  the  present  method  while  the  solid  curve  is  from  the 
reference.  The  results  agree  very  well. 
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C.  Buckling  of  plates  with  multiple  isotropic  stiffeners  and  sandwich 
panels  (RefsT  Iff,  20,  21 1  22) 

Figure  6.4  shows  seven  types  of  multiply  stiffened  plates  and  sandwich 
panels  that  are  used  to  correlate  with  the  present  method.  The  results 
of  correlations  are  given  in  Table  6.1  and  brief  descriptions  are  given 
below  according  to  its  order  in  the  table.  The  exactness  of  the  analysis 
in  these  references  varied,  but,  in  general,  the  translations  in  the 
Joints  are  neglected  while  moment  equilibriums  are  being  maintained. 

Panel  1  is  an  integrally  stiffened  plate.  When  the  length  of  the 
plate  is  so  proportioned  to  its  width  that  the  critical  mode  has  several 
axial  waves  and  the  buckling  load  is  not  affected  by  the  number  of  stiff¬ 
eners  which,  for  the  present  geometry,  is  about  four  or  five  stiffeners. 

Such  a  mode  is  called  local  mode  in  the  reference*^  and  the  stiffeners 
are  treated  as  discrete  plates.  The  result  shown  in  the  table  is  based 
on  an  infinitely  wide  plate  in  the  reference  and  for  a  six-stiffener 
plate  in  the  present  method.  When  the  panel  is  long,  such  as  a  25-inch 
long,  six-stiffener  plate,  the  present  method  shows  a  single  half-wave 
mode  and  the  buckling  load  is  22300  lbs.  where  stiffeners  are  treated 
as  discrete  beams.  It  differs  by  5%  with  the  result  from  Eq.  (A3) 
of  Ref.  20  which  treats  the  stiffeners  as  beams  with  their  elastic  proper¬ 
ties  smeared  onto  the  spacings. 

Panels  2  and  3  are  two  plates  with  Tee  and  L-section  stiffeners 
respectively.  They  are  from  the  same  reference  as  Panel  1  (Ref.  19) 
and  used  the  same  method.  Tee  and  L  sections  are  treated  as  discrete 
plate  elements  joined  togehter  with  translations  ignored.  In  the  reference, 
the  panel  is  taken  as  infinitely  wide  while  the  value  given  in  the  present 
method  is  for  a  six-stiffener  plate. 

Panel  h  is  an  isotropic  truss-core  sandwich  panel.  In  Ref.  21,  such 
a  panel  is  treated  by  assuming  various  buckling  modes  and  matching  moment 
equilibrium  in  joints  whose  translations  are  ignored.  The  lowest  load 
obtained  by  different  modes  is  taken  as  the  buckling  load.  This  approach, 
when  correctly  applied,  yields  very  good  agreement  with  the  present  more 
exact  method  as  can  be  seen  in  Table  6.1  where  two  cases  are  studied. 

The  only  difference  in  the  two  cases  is  the  thickness  of  the  wall  of 
the  core  (see  Fig.  6.1+) .  Case  A  is  for  a  core  with  equal  thickness  of 
core  web  and  the  face  sheet,  known  as  a  sandwich  where  core  restrains 
face  because  with  that  thickness  proportion,  the  face  sheet  relies  on 
core  webs  for  instability  support.  On  the  other  hand,  case  B,  where 
the  core  web  is  only  one-half  as  thick  as  the  face  sheet,  the  opposite 
is  true.  Figure  6.5  shows  the  mode  shapes  of  the  two  cases  by  calculating 
the  eigenvectors  of  the  buckling  analysis.  It  clearly  shows  that,  in 
case  A,  where  core-restrains-face,  the  face  sheets  near  the  two  free 
edges  of  the  panel  buckled  while  in  case  B,  the  face  sheets  are  thick 
enough  to  sustain  the  load  at  the  free  edges  of  the  panel  and  it  is 
the  core  web  near  the  center  of  the  panel  that  buckled  first. 
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Panels  5.  6  and  7  are  a  corrugated  sandwich  panel,  Zee  and  hat- 
section  stiffened  plates  respectively,  taken  from  Ref.  22.  The  method 
used  is  the  same  as  in  Ref.  19  where  the  stiffeners  and  the  web  of  the 
corrugation  are  treated  as  discrete  plates  Joined  together  and  displace¬ 
ment  of  the  Joint  is  neglected.  The  values  given  in  the  reference  are 
for  infinitely  wide  panels  while  those  of  the  present  method  are  fo~ 
finite  panels.  The  higher  buckling  loads  predicted  by  the  present 
method  are  due  to  the  fact  that  the  flanges  of  the  Zee  and  hat— section 
stiffeners  are  considered  as  a  part  of  the  plate  in  the  present  method 
while  in  the  references  the  bending  stiffness  of  the  plate  and  the 
flange  are  calculated  separately  and  then  added  together.  The  idealiza¬ 
tion  for  the  stiffeners  at  the  joints  as  treated  in  Ref.  22  is  closer 
to  a  riveted  attachment  while  that  of  the  present  method  is  similar  to 
a  bonded  connection.  It  is  interesting  to  note  that  tests  conducted 
in  Ref.  23  showed  that  the  buckling  stress  of  a  bonded  Zee-section 
stiffened  plate  is  19#  higher  than  a  corresponding  riveted  Z-stiffener 
plate.  Incidently,  this  value  is  quite  close  to  the  l8#  as  indicated 
by  Panel  6  of  Table  6.1. 


Table  6.1  Correlations  of  Analytical  Results  in  Literature  and  the 
Present  Method  for  the  Seven  Panels  Shown  in  Fig.  6.4 


Comparison  of  Buckling  Load 

Panel 

Remark 

Existing  * 
Reference 

Present 

Method 

Ratio 

1 

Ref.  19,  stiffener  as 
discrete  plate 

ks  =  1.87 

1. 86(m=6) 

1.01 

Ref.  20,  stiffener  as 
smeared  beam 

Pcr  =  23,1*50(1) 

22,432(1) 

1.04 

2 

Ref.  19,  local  mode, 
stiffener  as  plate 

kg  -  4.30 

4.25(5) 

1.01 

3 

Ref.  19,  Case  A, 

acr  =  21,400 

21,300(9) 

1.01 

Case  B, 

acr  =  42,000 

41,600(6) 

1.01 

Case  C, 

ocr  =  41,100 

41,900(7) 

1.00 

* 

Ref.  21,  Case  A,  core- 
restrains-face  type 

°cr  =  16 .920 

16,954(7) 

1.00 

Ref.  21,  Case  B,  face- 
restrains-core  type 

acr  =  6,070 

6,019(9) 

1.01 

5 

Ref.  22 

acr  =  20,400(5) 

19,550(5) 

1.03 

6 

Ref.  22 

acr  =  40,900(7) 

49,800(12) 

0.82 

7 

. 

Ref.  22 

aor  =  47,700(10) 

52,500(17) 

0.91 

*  k  is  a  buckling  load  parameter  which  is  proportional  to  the  buckling 
lSad,  PCf.  Is  the  total  load,  lbs.,  and  aCf.  is  the  stress,  pat. 
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From  the  previous  comparisons  with  existing  analytical  solutions 
for  isotropic  or  orthotropically  reinforced  structural  sections  and 
stiffened  plates,  one  can  see  that  the  correlations  are  usually  within 
a  few  percent.  In  a  few  cases  where  discrepancies  are  not  negligible, 
possible  causes  for  such  deviations  are  explained.  It  seems  that  the 
existing  approximate  methods  are  quite  accurate  and  in  most  cases  should 
require  less  computation  time  in  execution  than  the  present  method. 
However,  some  degrees  of  uncertainty  are  always  present  in  such  approxi¬ 
mate  methods  which  might  be  due  to  mode— guessing,  less  exactness  in 
matching  elements,  omission  of  coupling  between  in-plane  displacements 
and  bending  loads,  etc.  In  such  instances,  a  more  rigorous  theory 
like  the  present  method  should  be  employed. 

7.  TEST  CORRELATIONS  WITH  BORON  REINFORCED  PLATES, 

SECTIONS  AND  STIFFENED  PLATES 

Correlations  with  test  data  from  literature 

A.  Buckling  of  axially  loaded  all-composite  laminated  flat  plates 
(Refs.  24  and  25) 

Four  specimens  which  have  simply  supported  loaded  edges  are  taken 
from  Ref.  2k  to  correlate  with  the  present  theory.  The  physical  proper¬ 
ties  of  boron  tape  are  given  in  Ref.  2h.  The  test  specimens  are  20- ply 
plates  and  the  nominal  thickness  of  boron  tape  is  0.0053  inch. 

Table  7*1  shows  test  correlations  with  the  present  analysis  for  four 
plates  of  size  11  x  9.95  in.  with  the  length  11  in.  along  the  loading 
direction.  Correlations  are  reasonably  good. 


Table  7.1  Test  Correlations  for  All-Boron  Composite  Plates  (Ref.  2k) 


Condition 

Fiber 

Buckling  Load 

,  Nx,  lb/in. 

of  Unloaded 

Direction 

Analysis 

Sides 

ICO  4 

Other 

I  I 

Theories 

Value 

Rat  io( test/ analysis) 

1*04 

free 

parallel  to 
load  (0°) 

199 

215* 

216** 

206.5 

0.97 

405 

free 

perpendicular 
to  load  (90°) 

23.3 

23.7* 

23.6** 

22.6 

1.01 

4o4 

simple 

support 

parallel  to 
load  (0°) 

271 

292* 

299** 

285*** 

286.5 

0.95 

405 

simple 

support 

perpendicular 
to  load  (90°) 

251 

223* 

226** 

210*** 

217.0 

1.16 

*  Orthotropic  theory 
**  Rayleigh-Rit z  method 
**#Galerkin  method 
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Good  correlation  has  alao  been  obtained  from  Ashton  and  Love's  test 
(Ref.  25).  The  boron  fiber  orientation  (0°  and  combinations  of  0  and 
90°)*of  the  specimens  and  the  buckling  load  are  shpwn  in  Table  7.2. 

The  elastic  constants  of  boron  and  the  measured  thicknesses  of  plates 
that  were  used  in  calculations  are  given  in  Ref.  25.  The  results  showed 
that  the  maximum  error  is  less  than  3%  for  the  five  cases  where  load 
is  applied  at  0°  with  respect  to  fiber  orientation  and  much  larger 
discrepancies  for  the  other  five  cases  where  load  is  applied  at  90 
to  the  fiber.  Test  results  were  obtained  by  the  Southwell  plot.  The 
90°  tests  used  the  same  plate  after  the  0°  test  was  finished  by  turning 
the  same  plate  around  90°.  The  inferior  correlations _ in  the  cases  of 
90°  specimens  may  have  been  caused,  as  was  suggested  in  Ref.  25,  by 
effects  due  to  boundary  constraints.  The  test  specimens  are  20-ply 
square  plate  measured  10-inches  across  between  edge  supports.  The 
loaded  edges  are  clamped  and  unloaded  edges  are  simply  supp 


Table  7.2  Test  Correlations  for  All-Boron  Composite  Plates  (Ref.  25) 


Load  Applied 
Parallel  to  0°* 


Fiber 

Direction 

Test 
(lbs. ) 

Present 

Analysis 

Ratio 

0°  -  90° 
alternate 

11,300 

11,700 

0.97 

same 

9,500 

9,170 

1.03 

same 

7,400 

7,300 

1.01 

all  0° 

12,400 

12,100 

1.02 

all  0° 

13,700 

13,350 

1.03 

Load  Applied 
Perpendicular  to  0° 

Test  Present  Ratio 

(lbs.)  Analysis 

9.600  9,900  0.97 


7,700  1.17 
6,170  1.16 
3,040  1.43 
3,360  1.25 


#0°  refers  to  the  direction  of  fibers. 

Correlations  with  Boeing’s  tests  ofbucklinfi,  of  axially  compressed  composite 
reinforced  plates,  structural  sections  and  stiffened  pla  es. 


The  physical  properties  used  in  the  analysis  are: 


toron-e-poxy  BP  90 


Eu  =  29.117  x  10b  psi 
E22  =  2.341  x  106  psi 
G  =  0.75  x  106  psi 

V12  =  °* 2k61 

3 

Density  =  .072  Ib/in. 


Titanium  6A1-^V 
E  =  16.4  x  10^  psi 

G  =  6.2  x  10^  psi 

v  =  0.3 


3 

Density  =  0.158  lb/in. 


(7.1) 
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For  boron  composite,  and  Eg2  are  calculated  from  a  program  based 
on  Ref.  26.  The  G  value  from  the  program  seemed  too  high  and  the  present 
value  of  0.75  x  10^  is  the  same  as  that  in  Ref.  25.  Both  E22  and  G 
used  here  are  higher  than  that  in  Eq.  (6.l)  which  are  from  Ref.  18. 

For  aluminum  alloy  (7075— T6),  E  =  10.5  x  10^  psi  is  used.  Elasticity 
properties  of  adhesive  material  which  bond  composite  laminates  to  the 
metal  base  are  ignored  but  its  thickness  is  retained  in  the  analysis. 

B.  Test  correlations  of  composite  reinforced  flat  plates  (Boeing) 

Test  correlations  have  been  made  with  recent  Boeing  tests  on  composite 
plates.  A  total  of  1+8  specimens  were  tested  in  which  half  is  symmetri¬ 
cally  laminated  and  the  other  half  is  unsymme trie ally  laminated.  The 
loaded  edges  are  clamped,  the  unloaded  edges  are  free  or  simply  supported. 
Measured  between  clamps  and  knife-edges,  the  plate  is  9.0  x  2.98  inches. 
Thickness  of  the  specimen  is  from  0.097  to  0.127  inches.  Details  of 
the  geometry  are  given  in  the  Appendix. 

Test  results  and  correlations  with  analysis  are  given  in  Table  7«3. 

For  specimens  with  unloaded  sides  free,  the  calculated  buckling  load 
which  is  for  models  with  simply  supported  loaded  edges  is  multiplied 
by  a  factor  of  four  to  get  the  corresponding  value  for  the  clamped- free 
case.  That  a  factor  of  four  is  theoretically  correct  to  be  used  in 
the  present  case  is  verified  numerically  by  making  an  analysis  of  a 
plate  with  loaded  edges  clamped  and  unloaded  sides  simply  supported 
and  then  increase  its  width  to  extremity  (numerically,  from  the  original 
2.98  in.  to  500.0  in.).  The  buckling  load  of  such  a  wide  plate  should 
be  almost  the  same  as  a  clamped— free  plate.  The  result  verified  that  a 
factor  of  four  to  convert  the  simple-free  plate  to  a  corresponding  clamped- 
free  plate  is  numerically  correct. 

In  Table  7,3,  errors  based,  on  two  methods  of  evaluation  are  given. 

Method  1  is  based  on  using  the  measured  length  between  the  clamped 
edges  which  is  9*0  inches  as  the  theoretical  length  of  the  plate  in  the 
analysis.  This  implies  that  the  test  rig  is  theoretically  perfect  in 
making  the  ends  clamped.  The  average  errors  for  the  clamped-simple 
specimens  are  12 . 7$  and  13.6*,  and  for  clamped-free  specimens  are 
19.7#  and  38/5.  Since  in  the  test  set-up,  the  ends  of  the  specimen  are 
clamped  between  two  thick  blocks  which  are  then  clamped  to  the  testing 
head,  it  is  possible  that  perfect  clamped  conditions  might  not  have 
been  achieved  in  the  tests.  If  this  assumption  is  true,  one  may  postu¬ 
late  that  an  end-fixation  factor  of  four  used  for  the  clamped-free 
specimens  should  be  revised.  Based  on  the  24  test  results  of  the 
clamped-free  specimens,  an  average  end- fixation  factor  is  found  to  be 
3.14.  Assuming  that  the  buckling  load  is  proportional  to  the  end- 
fixation  factor  and  inversely  proportional  to  the  square  of  the  length 
of  the  plate  (true  for  the  classical  isotropic  plate  theory) ,  a  modified 
length  of  10.15  inches  is  found  (10.15  =  9  x  (4/3. 14)1'2) .  Based  on 
this  modified  length,  the  errors  of  correlations  are  given  in  Method  2 
of  Table  7.3.  The  average  errors  for  the  clamped-simple  specimens  are 
reduced  to  9%  and  5.5$;  and  for  the  clamped-free  specimens,  they  become 

9.9%  and  13.5*. 
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Table  7*3  Test  Correlations  of  Symmetrically  and 
Unsymme tri cally  Laminated  Composite  Plates  (Boeing) 


Test 

Specimen 

Boundary  Conditions : 

Loaded  Edges  Clamped , 

Others  Simply  Supported 

Boundary  Conditions: 

Loaded  Edges  Clamped, 

Others  Free 

Test,  Nx 
(lbs/in. ) 

Analysis 

Test,  N 

*  X 

(ibs/in. ) 

Analysis 

Method  1# 
Error 

w 

Method  2 ♦♦ 
Error 
<*) 

Method  !♦ 
Error 
(%) 

Method  2»* 
Error 
(%) 

nd 

<v 

8A-1 

4.0 

675 

-22.1 

4.0 

n 

-2 

4.6 

692 

-21.4 

4.6 

a 

-3 

3,020 

-11.7 

21.6 

733 

+  0.1 

21.6 

1 

8B-1 

3,690 

-14.6 

3.2 

751 

-23.0 

3.2 

s 

-2 

3,825 

-13.5 

0.6 

751 

-26.2 

0.6 

H 

d 

-3 

3,360 

-22.5 

5.5 

751 

-20.4 

5.5 

CD 

o 

•H 

8C-1 

3,690 

-  4.8 

7.0 

725 

-18.1 

7.0 

-P 

-2 

3,290 

-14.3 

7.7 

711 

-17.6 

7.7 

i 

-3 

2,750 

-36.6 

5.0 

691 

-21.0 

5.0 

w 

8D-1 

4,300 

-  5.6 

5.5 

835 

-20.2 

5.5 

a 

-2 

4,100 

-13.4 

4.1 

847 

-21.8 

4.1 

-3 

4,100 

-13.4 

2.4 

832 

-24.0 

_ 

2.4 

Average  error 

-12.7$ 

+9.07$ 

Average  error|l9.7#| 

***  +5.93$ 

8E-1 

3,830 

-10.5 

-  0.1 

1,340 

-13.4 

11.7 

<u 

-p 

-2 

4,025 

-12.6 

-  2.1 

l,l4o 

-43.2 

-12.6 

aJ 

d 

-3 

3,560 

-27.4 

-15.5 

1,010 

-61.6 

-27.0 

1 

8F-1 

5,280 

-  4.1 

7.1 

1,230 

-53.5 

-20.7 

-2 

3,840 

-32.8 

-18.3 

1,310 

-32.8 

-  4.3 

>> 

r*H 

-3 

4,070 

-22.9 

-  9.3 

i,4io 

-20.9 

5.0 

•d 

o 

8G-1 

5,175 

-  7.7 

3.4 

1,110 

-66.1 

-31.0 

*rl 

u 

-2 

4,740 

-13.6 

-  2.1 

1,110 

-60.0 

-25.7 

«P 

-3 

4,470 

-  6.3 

4.2 

i,i4o 

-36.1 

-  6.9 

i 

8h-i 

3,960 

-  9.2 

-  2.1 

940 

-26.8 

-  0.3 

w 

-2 

4,260 

-  7.6 

—  0.2 

970 

-32.4 

-  4.2 

-3 

4,025 

-  9.0 

-  1.8 

1,095 

-10.7 

12.7 

Average  error 

-13.6$ 

15.52$ 

Average  error  -38.1# 

|l3.5$ 

♦Method  1  used  the  measured  length  between  the  clamped  edges  which  is  9*0  inches 
as  the  length  of  the  plate  in  the  analysis.  Error  =  1  -  (Analytic  result /Test ) . 

♦♦Method  2  used  the  modified  length  which  is  calculated  as  10.15  inches  in  the 
analysis.  Error  -  1  -  (Analytic  result/Test ) . 

♦♦♦Average  of  the  absolute  values  of  the  errors. 


The  considerable  discrepancy  between  test  and  analysis  shown  in 
Table  7«3  may  be  caused  by  (l)  the  constraints  provided  by  the  test 
setup,  (2)  the  non-uniformity  of  specimens,  and  (3)  uncertainties  in 
determining  the  buckling  load  from  the  load- deflect ion  curves.  For 
item  (l),  the  use  of  effective  length  as  suggested  in  Method  2  of 
Table  7.3  seemed  to  have  improved  the  correlations.  However,  it  is 
not  a  satisfactory  answer.  The  non-uniformity  in  specimens  with  respect 
to  the  thickness  distribution,  as  indicated  in  item  (2),  may  cause  a 
symmetric  specimen  to  become  un symmetric.  The  effect  of  eccentricity 
on  the  buckling  load  will  be  less  severe  for  plates  with  sides  simply 
supported  than  those  with  sides  free  because  the  latter  has  no  support 
from  the  free  sides  to  resist  in-plane  bending  moments  produced  by 
the  eccentricity.  In  item  (3),  which  appeared  in  some  test  results, 
it  is  difficult  to  determine  the  buckling  load  from  the  longitudinal 
load-deflection  curve.  Two  examples  which  showed  such  undeterminable 
curves  are  given  in  the  Appendix  together  with  the  Southwell  plot  which 
used  the  lateral  load-de  fleet  ion  curve  to  determine  the  buckling  load. 

C.  Test  correlations  of  Boeing  specimens  of  composite  reinforced  struc¬ 
tural  sections  ( Fi g.  7.l) 

Figure  7*1  shows  the  geometries  of  some  structural  sections,  including 
angles.  Zee’s,  hat-sections  and  Tee-sections,  reinforced  with  uni-direc- 
tional  boron  strips  and  rods.  These  section  specimens,  as  well  as  the 
specimens  of  stiffened  panels  to  be  described  in  (7»  D* )  a re  machined 
flat  at  the  loaded  edges,  and  placed  in  the  testing  machine  without  any 
mechanical  fastening  between  the  specimen  and  the  testing  machine.  The 
unloaded  edges  of  the  section  are  free.  In  the  analysis,  the  lips  of 
Zee-sections  and  the  tips  of  Tee-sections  are  idealized  as  beams  tfhile 
all  the  other  reinforced  parts  are  treated  as  plate  elements.  The 
junction  of  the  machined  Tee-section  of  configurations  91-1  and  9J-1 
has  been  treated  as  the  nodal  line  of  three  plate  elements  without  the 
fillet  and  also  treated  as  a  beam  element  with  the  fillet  considered 
and  joined  by  three  plate  elements.  The  difference  in  the  results 
obtained  by  the  two  different  idealizations  was  found  to  be  small. 

The  test  results  and  correlations  with  analysis  are  given  in  Table  7«^* 
It  can  be  seen  that  all  correlations  are  reasonably  good  except  9A-1  and 
9J-1.  For  specimens  of  9A-1,  the  low  prediction  might  be  the  result  of 
idealization  which  had  neglected  the  corner  curvature  of  the  angle- sect ion. 
If  a  curved  plate  element  is  added  to  the  present  analysis ,  the  quality 
of  the  analysis  could  be  improved.  For  specimens  of  9J~1*  ihe  errors 
might  be  caused  by  difficulties  in  determining  the  buckling  load  from 
load-deflection  curves  of  test.  It  seems  that  for  composite  reinforced 
plates,  there  is  an  initial  period  of  adjustment  of  a  macroscopic  nature 
in  the  initial  loading  stage.  As  a  consequence,  in  some  specimens  at 
least, the  intial  part  of  the  load-deflection  curve  behaves  abnormally 
and  is  not  dependable.  If  a  specimen  is  such  that  its  buckling  load  is 
much  smaller  than  the  ultimate  load,  as  so  happened  in  the  case  of  9«J-1  » 
one  has  difficulty  in  reading  accurately  the  buckling  load  from  the 
test  curves. 
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Table  7.4  Test  Correlations  of  Boeing  Specimens  01 
Boron-Reinforced  Structural  Sections  (Fig.  7«1) 


r - - 

SPEC. 

NO. 

LENGTH 

INCH 

TEST 

TEST 7  ELASTIC  BUCKLINO 

ANALYSIS,  BUCKLING  PREDICTION 

CONFIGURATION 

ULTIMATE 
LOAD,  LBS . 

LOAD,  LBS. 

AVERAGE 

LOAD,  LBS. 

RATIO 

TEST/AN ALY . 

REMARK 

Ei*?±r=xi: xtZT'rm  ~  t. : 

REINFORCED 

TITANIUM 

ANGLE 

9A-1 

-3 

-4 

12. 0 

ft, 040 
8,660 
8,380 
8,420 

3,300  (0.0012)* 
3,400  (0.0011) 
3,600  (0.0011) 
3,700  (0.0015) 

3,500 

2,390  (0.00051) 

1 

1.46 

Low  analytic 
vnlue  le  due  in 
neglecting  the 
corner  engle  of 
apeclmen. 

Better  when 
apecleen  le 
longer,  . 

j 

rT 

-7 

20.0 

4,500 

6,800 

6,600 

2,200  (0.00065) 
2,000  (0.00050) 
2,200  (0.00040) 

2,130 

2,150  (0.00046) 

1 

. J 

0.99 

j  “SErNFORCK) 
TITANIUM 
!  ANGLE 

1  98-1 

f~-l 

-2 

-3 

1  ',.6 

11,240 

10,040 

9,140 

5,166“ 

(4,800)** 

6,000 

5,900 

7,010  (0.00098) 

1 

0.85 

l 

REINFORCED 

TITANIUM 

ZEE  9C-1 

-1 

-2 

-3 

13.6 

27,700 

27,500 

26,800 

20,000  (0.0028) 
17,000  (0.0018) 
22,000  (0.0024) 

19,700 

18,180  (0.0020) 

2 

1.08 

_ _ _ 

REINFORCED 

TITANIUM 

ZEE  9D-1 

-1 

-■* 

12.0 

11,240 

10,040 

9,140 

9,260 

5.800  (0.0017) 

4.800  (0.0011) 

5,850  (0.0015) 
5,600  (0.0017) 

5.510 

5.71*)  <0.00093) 

2 

0.96 

-5 

* 

20.0 

12,200 

11,500 

11,480 

— 

5,600  (0.0014) 
5,200  (0.0012) 
5,100  (0.0011) 

5.300 

5,660  (0.00092) 

3 

0.95 

REINFORCED 

TITANIUM 

HAT  9B-1 

-1 

-2 

-3 

13.6 

32,200 

30,500 

33,300 

32,200  (0.0037) 
30,500  (0.0030) 
33,300  (0.0033) 

33,300 

35,890  (0.0057) 

1 

0.95 

-1  end  -2 
delealnated  at 
the  given  loeda 

REINFORCED 

TITANIUM 

HAT  9r-l 

-1 

-2 

-3 

12.8 

6,900 

4,300  (0.0020) 

4,300 

4,590  (0.0016) 

5 

0.98 

-1  teated  to 

3,840  no  failure. 
-2  debond  at 

JiZOO  lbiu  r-^ 

REINFORCED 

ALUMINUM 

(FORMED) 

TEE  9G-1 

-1 

-2 

-3 

6.0 

20,300 

19,080 

22,320 

17,000  (0.0012) 
18,200  (0.0012) 
not  oertain 

17,620 

18,87 0  (0.0017) 

1 

0.93 

WSTNFQN0I5D 

ALUMINUM 

(FORMED) 

TEE  9H-1 

-1 

-2 

-3 

6.0 

20,940 

21,300 

22,320 

20,940  (0.0066) 
21,300  (0,0071) 
22,320  (0,0083) 

21,520 

23,520  (0.0041) 

1 

• 

0.92 

• 

REINFORCED 
TITANIUM 
(MACHINED) 
TEE  91-1 

-1 

-2 

-3 

-4 

10.0 

38,000 

37,200 

38,600 

37,000 

25,000  (0.0038) 
25,000  (0.0037) 
26,000  (0,0040) 
37,000  (0.0056) 

28,250 

28,830  (0.0041) 

1 

0.98 

REINFORCED 
TITANIUM 
(MACHINED) 
TEE  9J-1 

-1 

-2 

-3 

9.2 

11,040 

11,440 

10,260 

2,200  (0.00087) 
3,000  (0.00098) 
1,600  (0.00071) 

2,270 

1,850  (0.00067) 

1; 

1.23 

1  Axail  half-wave  nunber 


•  Unit  atrain,  in ./In.  ' 

**  The  original  tost  plot  of  load-va. -deflection  la  not  available  for  recheck  for  the  three  epeciaenB  of 
9B-1.  The  value  for  speclaen  -2  la  not  consistent  with  the  other  two  epeciaens  and  Is  not  included 
in  the  average* 


De  Tests  of  Boeing  specimens  of  composite  reinforced  stiffened  plates 

Figure  7.2  shovs  the  geometries  of  three  types  of  stiffened  plates, 
including  composite  reinforced  hat  and  angle-section  stiffened  plates 
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and  isotropic  honeycomb  sandwich  plates  reinforced  internally  with  boron 
composite  strips.  The  specimens  are  machined  flat  at  the  ends  and  placed 
in  the  testing  machine  without  any  fastenings.  The  unloaded  sides  of 
the  panel  and  stiffeners  are  free.  For  long  plates,  the  skin  is  instru¬ 
mented  for  elastic  buckling  and  the  skin  buckling  load  is  taken  for 
correlation  with  analysis.  For  short  plates,  only  ultimate  loads  were 
recorded  and  no  analytic  correlation  will  be  given.  In  the  analysis, 
the  deep  boron  reinforcement  of  the  stiffeners  in  plates  -A,  -B,  -0, 

-H,  -I  and  -J  is  treated  as  beams  connected  to  the  two  sides  of  angles 
as  plate  elements  while  in  -C,  -D,  -E  and  -F,  the  boron  strips  together 
with  the  immediate  skins  are  treated  as  plate  elements. 

From  Table  T«5,  it  can  be  seen  that  all  the  predicted  elastic  buckling 
loads  are  lower  than  the  test  ultimate  loads ,  except  specimens  -E  and  -F, 
which  are  honeycomb  sandwich  plates  whose  core  was  crushed  before  buckling 
occurs,  and  specimen  -I  where  the  two  numbers  are  extremely  close.  The 
polyimide  coat  seemed  greatly  strengthened  the  buckling  strength  of  the 
skin  in  specimen  -I.  In  general,  the  correlations  are  satisfactory. 

The  higher  predictions  are  due,  at  least  in  part,  to  the  fact  that  the 
skin  between  stiffeners  are  supported  laterally  only  by  frictions  developed 
during  the  test  while  in  the  theory  simple  support  of  the  skin  is  assumed. 


Table  7.5  Test  Correlations  of  Boeing  Specimens  of 
Boron-Reinforced  Stiffened  Plates  (Fig.  7.2) 


PART 

DESCRIPTION 

LOUTH 

TEST  ULTIMATE 

!  ANALYSIS,  BUCKLING  PREDICTION  I 

NO. 

INCH 

LOAD,  KIPS 

IBi 

TEST/ ANALYSIS 

U-A 

REINFORCED  ALUMINUM 
ANGLE,  LONG  PLATE 

33.7 

180  (0.0048)* 

129  (0.0028) 
akin  buckled 

150.8  (0.0033) 

12 

0.85 

11-B 

REINFORCED  ALUMINUM 
ANGLE,  SHORT  PLATE 

15.0 

210  (0.0056) 

N.  A.** 

154.1  (0.0033) 

6 

11  -C 

REINFORCED  TITANIUM 
HAT,  LONG  PLATE 

33.7 

155  (0.0052) 

skin  buckled  between 

25  and  50  kips.  Take 

30  kips  for  correlat¬ 
ion. 

33.5  (0.00088) 

14 

0.90 

11-0 

REINFORCED  TITANIUM 
HAT,  SHORT  PLATE 

15.0 

181  (0.0061) 

N,  A. 

33.8  (0.00089) 

6 

11 -E 

REINFORCED  TITANIUM 
HONEYCOMB  SANDWICH, 
LONG  PLATE 

33.7 

111  (0.0028) 

core  crushed 

150.4  (0.0038) 

1 

11 -r 

REINFORCED  TITANIUM 
HONEYCOMB  SANDilCH, 
SHORT  PLATE 

15.0 

211  (0.0056) 

eor#  failure 

200.5  (0.0051) 

1 

11-0 

REINFORCED  TITANIUM 
ANGLE,  LONG  PLATE 

33.7 

356  (0.0043) 

257  (0.0032) 
skin  buckled 

343.4  (0.0038) 

12 

0.75 

11-H 

REINFORCED  TITANIUM 
ANGLE,  SHORT  PLATE 

15.0 

461  (0.0056) 

N.  A. 

345.9  (0.0038) 

U-I 

j 

REINFORCED  TITANIUM 
ANGLE,  LONG  PLATE 
(POLYIMIDE  COATED) 

33.7 

350  (0.0042) 

325  (0.0026) 
skin  .buckled 

351.4  (0.0038) 

0.93 

11 -J 

REINFORCED  TITANIUM 
ANGLE,  SHORT  PLATE 
(POLYIMIDE  COATED) 

15.0 

232  (0.0028) 

1 

N.  A. 

349.9  (0.0039) 

^  Axial  half-war.  nunb.r 


•  Unit  strain,  in. /in. 

"  The  short  specimens  were  instrumented  for  crippling;  data  of  skin  buckling  was  not  available. 
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8.  DISCUSSIONS 


Some  points  of  interest  that  came  up  during  the  numerical  study  of  the 
method  will  be  discussed  in  the  following. 


Effect  of  layer  arrangement  of  composite  reinforced  plates 


A  simple  calculation  which  shows  the  effect  of  lamina  configuration 
on  relative  strength  and  weight  of  a  plate  under  axial  compression  is 
shown  in  Table  8.1. 


For  the  four  boundary  conditions  studied,  the  favorable  layer 
arrangement,  except  the  free-free  case,  is  to  place  composite  laminates 
between  metal  layers.  In  such  configuration,  weight  ratio  is  13% 
while  strength  ratio  is  95%  in  comparison  with  the  all-titanium  plate. 
When  the  composite  is  placed  outside,  strength  ratio  is  dropped  to 
about  50%.  A  possible  explanation  is  that  when  the  axially  compressed 
plate  is  buckled,  its  outer  elements  farther  away  from  the  mid-plane 
are  subjected  to  more  severe  extensional  and  twisting  strains  than 
the  elements  which  are  closer  to  the  mid-plane  and,  consequently, 
adequate  face  layers  should  have  high  G  and  E22  values  as  well  as 
value.  For  the  center  layer,  since  most  of  its  deformation  is 
relatively  of  the  axial  extensional  type,  its  E^  may  be  more  influ¬ 
ential  to  affect  the  strength  of  the  plate  than  its  E22  and  G.  This 
may  also  explain  the  low  strength  ratio  of  plate  No.  4  in  Table  8.1, 
since  in  that  case  the  composites  in  one  side  of  the  face  of  the  plate 
still  subjected  to  twisting  ■''n  buckling  deformation. 


Table  8.1  Uniaxial  Compressive  Buckling  of  Laminated  Plates 
(all  0°  fiber)  with  Various  Thickness  Ratios 
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B.  Effect  of  in-plane  constraints  along  the  unloaded  sides  of  a  stiffened 
plate 

The  effect  of  in-plane  displacements  along  the  unloaded  sides  of  a 
stiffened  plate  can  be  seen  from  a  calculation  made  for  the  single  stiffener 
plate  of  Fig.  6.2.  In  the  related  reference^®,  the  in— plane  constraints 
at  the  two  sides  are  =  Np2  ~  0  where  is  the  shear  force.  If  the 

condition  Wj.j?  =  0  is  replaced  by  u  -  0  which  means  the  displacement  in  the 
shear  force  direction  is  prevented  instead  of  the  shear  force  itself ,  the 
present  analysis  shows  that  the  buckling  load  will  be  increased  considerably. 
Table  8.2  gives  the  values. 


Table  8.2  Effects  of  In-Plane  Constraints  at  the 
Unloaded  Sides  and  Different  Stiffener  Idealization 
for  the  Stiffened  Plate  of  Fig.  6.2 


6  Value  in 
Fig.  6.2 

In-Plane  Boundary 
Constraints 

Stiffener 

Idealization 

Relative  Buckling 

Load 

N22  =  N12  =  0 

as  beam  (Ref.  18) 

1.00  (N/N  =  0.93)* 

P 

«22  *  »12  -  0 

as  beam 

1.00 

=  u  =  0 

as  beam 

( present 

1.24 

22 

analysis) 

W22  =  u  =  0 

as  plate 

1.54 

W22  =  N12  =  0 

as  beam  (Ref.  18) 

1.00  (N/N  =  2.86)* 
eq 

&22  =  11 12  =  0 

as  beam 

1.00 

N00  =  u  =  0 

as  beam 

(present 

1.23 

22 

analysis ) 

N22  -  U  =  0 

as  plate 

1.32 

*E  A  /Ebt  =1.5,  see  Fig.  6.2 

S  5 


It  can  be  seen  that  for  the  prescribed  geometry,  buckling  strength  is 
increased  more  than  20#  by  preventing  the  sides  from  moving  in  the  direction 
of  the  load.  For  different  geometry,  however,  the  effect  may  not  be  so 
noticeable. 

C.  Difference  due  to  idealization  of  stiffeners  and  flanges 

Stiffeners  and  flanges  of  a  section  or  a  stiffened  plate  can  be  idealized 
into  either  beam  elements  or  plate  elements  depending  on  the  geometry  of 
the  member  itself  and  the  geometry  of  the  structure.  However,  different 
ways  of  idealization  may  produce  quite  substantial  discrepancies  in  the 
analytical  results.  Table  8.2  above  shows  the  difference  in  loads  by 
treating  the  single  stiffener  in  Fig.  6.2  first  as  a  beam  and  then  as  a 
plate.  Two  similar  instances  are  shown  in  Table  8.3  taken  from  structures 
of  Figs.  6.1  and  6.1*. 
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Table  8.3  Effect  of  Different  Ways  of 
Idealization  of  Stiffeners  and  Flanges 


Example  Descriptions 

Stiffener  or  Flange 

Idealization 

Buckling  Load 

Web  with  orthotropic 
flange .  Long-wave  mode 
described  in  Sec.  (6*  A) 
and  Fig.  6.1 

flange  as  beam  (Ref.  18) 

k  =  1.56 
w 

flange  as  beam  (present  analysis) 

1.58 

flange  as  plate  (present 
analysis) 

5.90 

Plate  stiffened  by  six 
isotropic  integral 
stiffeners.  Long-wave 
mode  described  in 

3ect.  (6,  C)  and  Fig. 
6.4,  panel  1. 

stiffeners  as  beams  (Ref.  20) 

P  =  23,450  lbs 
cr 

stiffeners  as  beams  (present 
analysis) 

22,347 

stiffeners  as  plates  (present 
analysis) 

29,468 

Such  large  discrepancy  as  shown  in  the  first  example  of  Table  8.2  in 
treating  the  flange  as  beam  and  as  plate  warrants  further  study  of  the 

matter. 

D.  Buckling  mode  shape  plots 

In  Fig.  6.5,  the  mode  shapes  of  a  truss  core  sandwich  panel  corresponding 
to  two  different  core  web  thicknesses  have  been  shown.  Such  plots  yield 
valuable  physical  insight  into  the  instability  mechanisms  of  structures. 

An  additional  example  is  given  here  to  demonstrate  the  usefulness  of  such 
mode  shape  plot  in  identifying  weak  members  in  a  structure  from  a  buckling 
point  of  view.  Fig.  8.1. a  shows  the  buckling  mode  shape  for  panel  1  in 
Fig.  6.4,  with  three  stiffeners.  It  is  seen  that  the  stiffener  is  weaker 
than  the  skin  from  the  point  of  view  of  buckling.  In  the  succeeding  figures 
of  Fig.  8.1,  the  stiffener  thickness  is  gradually  increased  and  the  skin 
thickness  correspondingly  decreased  so  that  the  total  cross-sectional  area 
(and  hence  the  weight)  is  the  same  in  all  cases.  It  is  seen  that,  for  the 
same  weight,  the  panel  of  Fig.  8.1.c  carries  the  maximum  load,  a  gain  of 
5k%  over  the  panel  of  Fig.  8.1. a.  Similar  studies  can  be  made  on  more 
complicated  cross-sections.  The  usefulness  of  the  mode  shape  plot  in 
assisting  design  and  in  structural  optimization  is  evident. 

9.  CONCLUSIONS 

The  present  analysis  for  the  instability  of  composite  plates,  sections  and 
stiffened  plates  with  composite  reinforcement,  is  an  exact  theory  in  the  classi¬ 
cal  sense.  The  connection  between  plate  elements  and  between  plate  and  beam 
elements,  with  the  effect  of  offset  included  in  the  analysis,  and  the  con¬ 
siderations  for  the  conditions  at  all  the  unloaded  edges  of  the  structure  are 
rigorous  and  consistent  with  the  linear  plate  theory  and  elementary  beam  theory. 
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Reasonably  good  agreement  with  existing  analytical  and  test  data  is 
obtained  from  the  results  of  the  correlation  study.  Some  scattering  in 
test  correlations  indicates  the  degree  of  uncertainties  in  dealing  with 
composite  materials  in  theory  as  well  as  in  techniques  of  fabrication  of 
the  composite  reinforced  stiffened  structures.  If  the  composites  are  so 
ideal  such  that  no  shear  strain  exists  in  the  direction  of  the  thickness, 
all  layers  are  homogeneously  orthotropic,  the  material  stress-strain 
follows  linear  Hooke's  law  until  buckling  occurs,  no  residual  stress,  no 
sliding  between  layers  and  no  separation  occurs  between  laminates  and 
metal,  then  the  mathematical  model  stipulated  for  laminated  composites 
in  the  analysis  is  exact  in  reality,  and  reasonable  correlation  between 
theory  and  tests  might  be  expected.  Further  complication  arises  from  the 
uncertain  properties  of  the  adhesive  layer  which  bond  the  composite  lamina 
to  the  metal.  When  the  thickness  of  the  adhesive  layer  is  of  comparable 
order  of  magnitude  as  the  thickness  of  the  composites,  it  could  be  impor¬ 
tant  to  include  the  elasticity  of  the  adhesive  layer  and  the  effect  of  the 
inter-lamina  shear  into  the  analysis. 

Since  the  present  method  combines  local  instability,  which  involves  only 
some  of  the  elements,  as  well  as  general  instability,  which  involves  the 
whole  structure,  the  solution  of  the  eigenvector  for  the  particular  eigen¬ 
value  (buckling  load)  is  useful.  This  eigenvector  capability  has  been 
included  and  it  would  be  useful  in  an  optimization  analysis  where  the 
weak  members  which  buckled  first  could  be  detected  and  reinforced. 

The  present  method  could  be  used  in  conjunction  with  existing  optimize 
tion  and  direct— search  programs  to  establish  mini mum- weight  configuration 
for  panels  with  composite  reinforced  stiffeners. 
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Fig*  5*3  A  joint  where  three  plate  elements  meet 

a  beam  element 


(1)  (2)  (U)  (7)  (8)  (9)  (ID  (HO  ( 22 ) ( 23 )  (25)  (28)  (29) 


Fig.  ^.1*  A  hat-section  stiffened  plate  with  four 


hat-sections 


BUC 


Fig*  6.1  Buckling  of  web- flange  with  and  without  composite 
reinforcement  at  the  flange  (Ref.  18) 


-  analytical  results  from  Ref.  18 

©  results  from  present  method 

X  =  axial  half-wave  buckle  length 

_  2  2  — 

k  =  N  b  / tt  D  where  N  »■  load  per  inch  on  web 

w  X  w  W  X 

Dw  =  bending  stiffness  of  web 
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E  A  /Ebt 
s  s 

Fig.  6.2  Discretely  stiffened  isotropic  plate  buckled  under  axial 

compression  -  single  shallow  eccentric  stiffener  (Ref.  18) 

E  =  Young’s  modulus  of  stiffener;  E  =  same  for  plate 
s 

A  =  cross  sectional  area  of  stiffener 
s 

N  -  load  per  inch  on  plate,  including  stiffener  load 

W  =  buckling  load  per  inch  on  an  equivalent  plate  of  equal 
eci*  weight  with  constant  thickness 


*axial  half-wave  number 


Fig.  6.3  Discretely  stiffened  plate  -  single  eccentric 
orthotropic  deep  stiffener  (Ref.  18) 


#  number  of  axial  half-wave 

Np  -  load  per  inch  on  isotropic  plate 

D  =  bending  stiffness  of  plate 

Geometry  at  points  1  and  2  on  Fig.  6.3:  (unit:  inch) 
Point  d  t  3  n  a  =■  length  of  plate 


s|*.«J  bflLL-4  bs  ibfi  bw  | 

rA  bf2  %f2 


• 

• — 1 
< 

al loy . 

6  stiffeners 

a  = 

12.3 

t  =  0.089 

s 

b 

s 

=  2.05 

t  =  0.058 
w 

b 

=  i.oo 

a  =  length  of 

w 

plate 

Al. 

alloy. 

6  stiffeners 

a  = 

15.0 

t  =0.080  t  =t 
s  f  w 

b 

s 

=  3.0 

t  =0.056 

w 

b 

w 

=  1.8 

bf=0.5U 

Al. 

alloy . 

6  stiffeners 

a  t 

=t„  b  b„ 

w  f  v  f 

A 

iu.lt 

.o4  1.92  .576 

B 

14.4 

.04  0.96  .288 

C 

16.0 

.048  1.50  .500 

Al. 

alloy. 

13  6o°-cells 

a  “ 

6.0 

t  =  0.02  for 
c  . 

case  A 

b 

s 

=  1.0 

=  0.01  for 
case  B 

+■ 

vt 

=  0.02 

Al. 

alloy. 

7  cells 

a  = 

8.96 

»fl  =  2.56 

b 

c 

=  1.06 

*f2  -  i-06 

t 

c 

=  .036 

tf  =  .048 

Al. 

alloy. 

6  stiffeners 

a  = 

16.0 

bf  =  0.500 

b 

w 

=  1.436 

t  =  0.080 

s 

bg 

=  1.876 

t  =  t  =  0.048 
w  f 

Al. 

alloy . 

5  stiffeners 

a  = 

20.0 

b  =  1.5 
s 

bn 

=  0.5 

b  =  1.436 

w 

bf2 

=  1.2 

Fig.  6.1  Geometry  of  seven  types  of  isotropically  stiffened  plates 

and  sandwich  panels  for  analytical  correlations  (unit:  inch) 
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0.02  INCH 


Case  A  of  Panel  4  in  Fig.  6.4.  Core-restrains-face  type 


Case  B  of  Panel  4  in  Fig.  6.4.  Face-restrains-core  type 


Fig.  6.5  Buckling  shapes  of  a  truss -core  sandwich  panel 

Case  A:  Core-restrains-face  type 

(Core  web  thickness  =  0.2) 

Case  B:  Face-restrains-core  type 

(Core  web  thickness  =  0.01) 

Note:  The  relative  amplitude  of  each  element 
of  the  buckled  shape  is  drawn  in  scale 
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( 9  A-l )  8-ply,  Tol+O  TITS'  .03 
(9B-1)  12-ply,  .063  1.10  .12 


t  b  r 

(9C-1)  12-ply,  TO  1.10  .12 

(9D-1)  8-ply,  .0U0  1.16  .08 


ANGLE  ZEE 


— 


(9E-1)  12-ply,  7o53  P3  OF  .125 
(9F-1)  5-ply,  .025  0.50  0.50  .05 


Fig.  7.1 


Geometry  of  test  specimens  of  boron-reinforced  titanium  or 
aluminum  alloy  structural  sections* 

Materials  given  in  Table  7*^*  Ply  thickness  of  boron  tape 
*  0,0053  inch  (aOMlnaX)  Adhesive  thickness  =  0*018  approx. 


^  metal  ^///^ 


boron  composite  layers 
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Fig*  7*2  Geometry  of  test  specimens  for  composite  reinforced 

sti ffened  plates  and  stiffened  honeycomb-sandwich  panel 
(Unit:  inch) 

Thickness  of  each  boron  tape  layer  =  0.0053  in. 
Thickness  of  adhesive  layer  =  0.013  approx. 


59 


Configuration 

Buckling  type 

Pcr.  Its. 

m 

ecr,  in. /in 

A 

Stiffener 

"buckled 

27,51+0 

6 

0.00317 

B 

Stiffener 

buckled 

38,31+6 

5 

0.00420 

C 

Skin  and  Stiffeners 
buckled 

1+2,1+91 

5 

0.001+90 

D 

Skin  near  edge 
fliirmort  buckled 

35,285 

7 

0.00407 

Note:  All  configurations  have  the  same  cross  sectional  area 
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APPENDIX 


SUPPLEMENTARY  ANALYSIS, 

ANALYTICAL  CORRELATIONS  AND  GEOMETRIC  DETAILS 
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INTRODUCTION 

The  Appendix  of  the  report  contains  some  peripheral  equations  which  are 
not  essential  in  the  understanding  of  the  analysis  but  would  complement 
the  equations  in  the  main  text  in  working  out  numerical  problems  which 
contain  plate  elements  with  uncoupled  features  (B^  *  0).  Also  included 
in  the  Appendix  are  a  few  additional  analytical  correlations  and  a  detailed 
geometry  of  the  Boeing  specimens  of  composite  reinforced  plates. 

1.  EQUATIONS  TO  IMPLEMENT  EQUATIONS  IN  THE  MAIN  TEXT 

1.1  Equations  for  Rectangular  Composite  Plates 

In  the  main  text,  equations  for  plates  with  elastically  restrained  sides 
are  given  in  Eqs.  (5-1)  to  (5.10a).  For  other  boundary  conditions,  corres¬ 
ponding  equations  are  given  below. 


Plates  with  all  four  sides  simply  supported 

Assume  that  the  origin  of  the  axes  is  placed  at  a  corner  of  the  plate. 


The  boundary  conditions  are: 

x  =  0  or  a 

y  -  0  or  b 

(i)  w  =  0 

(ii)  M-q  =  0 

(ili)  -  o  (“-1’ 

(iv)  v  =  0 

(i)  w  =  0 

(ii)  M22  “  0 

(iii)  u  =  0 

(iv)  N22  =  0 

(II-2) 

«  ^  • 

Displacement  functions  assumed  which  automatically  satisfy  the  above  condi 


tions  are: 


V 

=  W 
mn 

sinot 

sin3 

V 

=  Vmn 

sina 

cos  3 

u 

Umn 

cos  a 

sin  3 

where 

a=  mirx/a  and  f$  =  nny/b  and  m  and  n  are  half-wave  numbers. 

On  substituting  the  displacement  functions  in  the  equilibrium  equations, 
one  arrives  at  an  equation  similar  to  Eq.  (5*6)  from  which  the  buckling 
load  3n  is  obtained  explicitly  as  a  function  of  m  and  n.  Vary  m  and  n, 
the  lowest  "S^.1  thus  obtained  is  the  buckling  load. 

Plates  with  loaded  edges  clamped,  other  two  sides  simply  suppo.rtgci 

Assume  a  coordinate  system  such  that  the  x-axis  coincides  with  one  si  e 
of  the  plate  and  the  y-axis  is  at  the  middle.  The  boundary  conditions  are: 

_ +£  y  =  0  or  b 

A  «  2  T  "  “  ' 

(i)  w  =  0  (II-1*)  (i)  w  =  0  (II-5) 
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(ii)  w  x  *  0 

(iii)  MU  =  0 

(iv)  v  =  0 


(ii)  Mg2  =  0 

(II-4)  (con't)  (iii)  u  =  0 
(iv)  N22  =  0 


( II— 5 )  (con’t) 


Displacement  functions  assumed  are  (upper  rows  correspond  to  symmetric  mode 
and  lower  rows  to  antisymmetric  mode) : 


r  TT  /coshou  .  Q 
w  =  )  W.  (  .  ,  )  sing 

.  S  i  sinhot 
i=l 


a  =  Tip^x/a 


r  „  /Coshou  „ 
}  V.  (  .  ,  )  cosg 

f  i  sinhtt' 


g  =  n*ny/b 


(II-6) 


v  TT  /Sinhcu  .  0 
)  U.  (  ,  )  sing 

.  i  cosha 
i=l 


which  satisfies  the  boundary  conditions  at  y  =  0  and  b. 

Since  the  displacements  satisfy  the  conditions  of  symmetry  or  antisymmetry 
at  the  middle  of  the  plate  (x  =  0),  only  the  four  conditions  at  x  =  a/2  need 
to  he  considered  and  only  four  roots  of  Eq.  (5.6)  will  he  used.  The  other 
four,  say,  the  negative  roots,  will  he  ignored. 

The  UxU  buckling  determinant  (Eq.  (5.10))  for  the  symmetric  axial  mode, 
can  be  written  as 


cosh(  ttp^/2) 
pisinh(  7rp^/2) 


( II— 7 ) 


tAllL2i(pi/a)"A12Lli(ll/b)"Bll(pi/a)2+B12(n/b)^lcOSh(  ‘  • 

L^cosh(  irp^/2)  ... 


i-1 »2  ,3 


where  Ln  and  bgj.  are  given  in  Eqs .  (L;.8).  For  the  antisymmetric  mode,  one 
simply  changes  cosh  to  sinh  and  sinh  to  cosh  in  the  above  equation. 

When  the  coupling  matrix  Bij  is  zero,  it  is  possible  to  reduce  the  buckling 
determinant  to  2x2  since  only  the  first  two  boundary  conditions,  (i)  and  (ii) 
of  Eq.  (II-4)  need  to  be  enforced. 

Plate  with  loaded  edges  simply  supported,  one  side  simply  supported,  the 

other  side  free  .  t 

The  coordinate  axes  are  the  same  as  the  elastically  restrained  plate.  The 

boundary  conditions  are: 
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x  =  0  or  a 


y  -  -b/2 


(l)  w  =  0 

(ii)  Mn  =  0 

(iii)  Nn  =  0 

(iv)  v  =  0 


( II— 8) 


(i)  w  =  0 

(ii)  M22  =  0 

(iii)  u  =  0 

(iv)  Nqo  =  0 


(i)  M22  =  0 


(ii) 


+  2 


(iii)  N12  =  0 

(iv)  N22  =  0 


(II-9) 


(11-10) 


Displacement  functions  assumed  are  the  same  as  Eq.  (5.1).  Proceeding 
exactly  as  before,  the  resulting  8x8  determinant  can  be  written  as 


exp( 

1Wi/2)G2i 

exp( 

wpi  /  2 )  Gii 

exp( 

^i^^Ui 

exp( 

exp( 

-np./2) 

exp( 

-7ip./2)G5i 

exp( 

-npi/2)L2i 

exp( 

-irpi/2)GUi 

=  0 


i  =  1,2, ...8 


(11-11) 


where  ,  L2^  are  given  in  Eq.  (5.8)  >  Gii  to  Gl^  are  given  in  Eq. 
and  05^  is  defined  as 

G5i  =  -B12(m/a)L2i+B22Lli(pi/b)+D12(m/a)c;-D22(pi/b) 


(5.10a) 

(11-12) 


When  Bji  *  0,  the  order  of  the  buckling  determinant  is  reduced  to  Hx4, 
since  the  boundary  conditions  that  need  to  be  enforced  are: 


y  »  -W2 


y  °  +E/& 


(1)  w  =  0 

(2)  M22  =  0 


(1)  Mg2  »  0 

3^22  ^12 

<2>  -#  *  2  TT-  ‘  0 


(II-I3) 
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Plate  vith  loaded  edges  simply  supported,  the  other  two  sides  free 
The  loaded  plate  has  the  same  coordinate  system  as  the  elastically 
restrained  plate.  The  boundary  conditions  are: 

x  -  Q  and  a  y  =  tb/2 


(i )  v  =  0 

(i)  m22  =  0 

(li)  Mu  =  0 

,  ...  m22  „  8M12  . 

(ll)  3 y  *  2  3*  ■  0 

(iii)  Ni;l  =  0 

(II-14) 

(iii)  »12  =  0 

(iv)  v  =  0 


(iv)  N22  =  0 


Displacement  functions  assumed  are  (the  upper  rows  correspond  to  symmetric 
mode  and  the  lower  rows  to  antisymmetric  mode): 


v-  y  W.(C03^)  Sine 
i  sinhe 

4 

v  =  l  vi( 

♦  _ _  *1 


i=l 


sinhfl 

cosh$ 


)  sin  a 


a  =  mirx/a 


6  =  "P^y/b 


< 11-16) 


r  TT  /COShftx 

u  *  )  U.  4  cos  a 

i  sinhft 

The  buckling  determinant  is  formed  from  the  enforced  boundary  conditions 
(Eq.  (11-15))  using  the  four  positive  roots  from  Eq.  (5.6).  The  resulting 
4x4  buckling  determinant  is  (for  the  symmetric  mode) 


cosh(  •npi/2)G5i 
sinh(irpi/2)Gli 
cosh(iTp^/2)  °4i 
sinh(irpi/2)G3i 


=  0 


i=l,2.3,4 


(II-17) 


where  ,  G3i,  are  given  in  Eq.  (5.10a)  and  G5i  is  given  in  Eq.  ( 11-12). 

When  B^j  =  0,  the  buckling  determinant  is  reduced  in  size  to  2x2  since 
only  the  first  two  boundary  conditions  of  Eq.  ( 11-15)  need  to  be  considered. 


1.2  Equations  t'or  structural  Sections  and  Stiffened  Plates 

Equations  necessary  for  the  analysis  involving  plate  elements  with  -  0 
are  complete  in  the  main  text.  The  corresponding  equations  for  the  simpler, 
uncoupled  cauen  of  By  =  0  will  be  given  here  for  reference. 
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Flat  plate  element-forces  and  displacements  in  local  coordinates 
For  the  case  of  Bij  —  Q*  the  displacement  functions  assumed  are. 


w  =  I  W 
i=l 


sina 


-  I  V- 


i=l 


u  =  I  U  e 
i=l 


g. 


sina 


coso 


(11-18; 


where  a(  a  =  nnrx/a)  is  the  same  as  before.  The  symbol  g  used  in  Eq.  (5.18) 
is  changed  to  gp  and  &2  or(^er  to  distinguish  the  characteristic  roots 
obtained  in  bending  and  in  stretching  which  are  now  uncoupled.  They  are 
defined  as : 


P1  =  Puiiry/a 
B2  =  Pwi^/a 


(11-19) 


'The  displacement  functions  after  satisfying  the  equilibrium  equations 
(through  Eqs.  (4.5)  and  (4.6))  reduce  to  the  determinantal  equation  of 
Eq.  (5.4)  where  R13,  R23,  R31  and  R32  are  now  zero.  This  results  in  two 
separate  equations.  The  first  one  is 


*11 


21 


R 


12 


22 


=  0 


(11-20) 


These  four  elements  are  the  same  as  in  Eq.  (5*5)  with  a  replacing  b  and 
Pui  replacing  p3  which,  after  expansion,  becomes 

K,p\  *Kn*o 

U4  ill  u2  U1  uu 

which  yields  four  roots  of  pu^  to  he  used  in  u  and  v  displacements  in 
Eq.  (11-18) .  Vi  in  Eq.  (Il-lo)  can  now  be  expressed  in  terms  of  Ui  ass 

(11-22) 


V.  =  L-.U. 
1  3i  1 


where 


L3i  =  “*11^*12  =  "R21^R22 
The  second  equation  is 

*33  =  ° 


(11-23) 


( 11-24) 


where  E™  is  the  same  as  in  Eq.  (5.5)  vith  a  replacing  b  and  Pvi  replacing 
.  It  is  a  fourth  order  polynomial  in  even  powers  of  Pvri-  The  four  values 
of  pwi  from  Eq.  (11-24)  are  to  be  used  with  w  displacements  in  Eq.  (U-lOi. 
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Using  these  p^  and  p^  values,  the  three  displacements  in  Eq.  (5-13^ 
can  he  expressed  in  terms  of  and  as: 


r  ?  r 

w  »  i  e  '  sina  -  l  uj^  W  sina 


4  61  4 

u  =  7  U.  e  cosa  =  T  n  .  U.  cosa 
i=l  1  i=l  ^  1 


(11-25) 


P1 

=  y  U.  .  e  sina  =  Y  p  .  U.  sina 
>.131  >,  ui  i 

i=l  i=l 


and  by  differentiation. 


p2 

wfy  *  I  Wi(npwi/a)e  sina  =  £  4>wi  W.  sina 


(II-26) 


Using  the  above  displacements,  the  four  stress  resultants  in  Eq.  (5-12) 
can  he  written  as  (through  Eqs.  (4.5)  and  (4.6)): 

Qg2  3  (i»12(m/a)2(pwi/a)-D22lpwi/a)3+UD66(m/a)f-(pwi/a)}Wi7i3e  2sina 


=  y  q^^W^sina 
i=l 


( 11-27 ^ 


k  9  p  p  Bp  It 

Ipp  3  l  {D12(m/a)  -D22(pvi/a)  ^ir  e  sina  =  £  (m22)wiWisinc 

^  i=l  i=l 

(11-28) 


k  gl  r 

N12  =  l  {Agg[(pui/a)+(m/a)L3i}Uiire  cosa  =  £  (n^J^^cosa 

i=1  i=1  (II-29) 

4  h  r 

«22  «  l  {-A12(m/a)+A22(pi,4/a)L3i}Uiire  sina  =  J  ^n22^uiUi8in° 

(II-30) 

Transformations  to  the  reference  plane 

When  Bij  =  0 ,  corresponding  to  Eqs.  (5.23)*  one  has 


w  =  y  (u>  .  +  y  4  .  )W.  sma 
s  L  '  wi  o  wi  i 
i=l 


(II-31) 
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(11-32) 


4 

w  =  y  <i>  •  .  W.  sina 
s,y  i^1  Ywi  i 


u  =  y  [n  .U,-z  (niTT/a)cj  .W.+y  (up  ./a)n  .U.Jcosa 
s  .4  ui  i  o  wi  i  o  *ui  ut  l 


i=l 


v  =  y  (p  .U.-z  d>  .W.)sina 
b  Mui  i  oYwi  1 


Or  v/ritten  in  matrix  form 


ws  ' 


s.y 


U 


f  x*  1 

I  (4x8)  I 


i.e. 


{dpg}  =  (Xj]  (R*) 


{R*}  is  self-evident 


W, 


U. 


u, 


The  boundary  forces,  for  B^j  =  0,  are 

4 

(q22}s  =  .n,wvvVui(m*/a)ui]sina 


(M22}s  =  .4[(m22)wiWi+yoWZo(n22)uiUi]sina 
1=1 

4 

<»12>,  =  j1(n12>ulUi“s'' 


4 

(N22)s  =  l  ^n22)ui+yo(ni2)ui(mir/a)]UiSina 
i=l 

Or  written  in  matrix  form: 


* 

1  r  1 

WM 

*xw 

I 

K 

<"l2>s 

OO 

X 

J- 

II 

V1 

1 

>  4  - 

U4 

(11-33) 

(11-34) 

(11-35) 

(11-36) 

(11-37) 

(11-38) 

(11-39) 

(II-40) 

(II-Ulj 
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i.e. 


(11-42) 


{fpsJ  =  [X*]  {R*} 

Flat-plate  element  -  forces  and  displacements  in  global  coordinates 
~  For  Bj,  =  0,  equations  similar  to  Eqs .  (5.33),  (5.3M»  ( 5  .  3^ ) and. 

(5.37)  are 

(dJG)  =  [Td]  [xj+]  {R*}  =  [X*+]  {R^}  y  =  +b/2  (11-43) 

{<^G)  =  [Td]  [X*-]  {R*}  =  [X*-]  {R*}  y  =  -b/2  ( 11-44) 

(fJG)  =  [Tf]  [X*+]  {Rj}  =  [Xj+]  {R*}  y  =  +b/2  (11-45) 

{fPG}  =  [Tf]  [X2~]  Ul}  =  [X4~]  {R1}  y  =  “b/2  (11-46) 

where  the  matrix  [x|+]  is  the  matrix  [X*]  with  its  y  replaced  by  +b/2,  etc. 

The  equations  for  a  plate  element  in  global  coordinates  are  complete  with 
the  addition  of  Eqs.  (11-43  to  46)  which  will  be  used  for  a  plate  element 
when  its  =  0.  In  such  case,  one  simply  uses  {R*},  [X^+],  ^"]t  [X*+] 

and  [XjJ“]  to  replace  (R^ ,  [Xjj],  [X~] ,  821(1  respectively. 

The  boundary  conditions  corresponding  to  a  plate  element  with  B-j^  =  0 
can  be  easily  derived  from  the  above  equations  and  will  not  be  presented 
here, 

2.  ADDITIONAL  ANALYTICAL  CORRELATIONS 

Two  more  correlations  of  the  present  method  with  analytical  results  from 
existing  literature  will  be  given  here. 

2.1  Buckling  of  an  Angle  Section  (Ref.  II-l) 

Figure  11,1  shows  the  buckling  load  of  a  column,  its  cross  section  is  an 
angle.  Curves  A  and  B  are  results  taken  from  Ref.  II-l.  The  former  is 
obtained  by  assuming  that  the  angle  is  made  of  two  plates,  each  one  is  simply 
supported  at  three  sides  and  the  fourth  side  is  free.  Correlation  with  the 
present  method  is  good  even  though  the  three-sides-siraply-supported  assump¬ 
tion  does  not  satisfy  exactly  the  continuity  conditions  along  the  common 
boundary  of  the  two  plate  elements.  The  second  curve, B,  which  treats  the 
column  as  an  Euler  column  with  pin  ends,  is  almost  one-half  the  value  obtained 
by  the  present  method  which  assumes  a  line  support  at  the  ends  with  simple 
support  conditions  rather  than  a  pin-end  phenomenon.  Considering  that  an 
Euler  column  with  clamped  ends  has  a  buckling  load  four  times  the  value  of 
its  pin-ended  counterpart,  the  result  from  the  present  method  which  is  about 
twice  the  pin-ended  column  solution  seems  reasonable. 

2.2  Buckling  of  Beaded  or  Lipped  Isotropic  Sections  and  Plates  (Refs.  9* 

I 1-2,  19,  22) 

Geometries  of  some  isotropic  sections  and  plates  with  beads,  flanges  and 
lips  are  shown  in  Fig.  II-2.  These  structures  and  the  corresponding  analytical 
results  of  buckling  are  taken  from  Refs.  9 9  II-2,  19  and  22  where  beads  and 
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lips  are  treated  as  beams.  The  major  difference  between  the  present  method 
and  these  references  is  the  omission  of  continuities  in  in-plane  displace¬ 
ments  at  common  boundaries  and  neglect  offset  between  plate  and  beam  elements 
in  these  references.  Table  II. 1  shows  the  results  of  correlations  which  seem 
to  be  reasonably  good. 


Table  II. 1  Correlation  of  beaded  or  lipped 
isotropic  sections  and  plates 


Buckling  Load,  ^j.1*  i 

on  the  plate 

element 

without  including  load  on  stiffener 

Section 

Dnowii 

figure 

From  reference 

Present 

Ratio 

II-2 

Buckling  load 

Ref.  No. 

Method 

Beaded  plate 

-a 

12,280  ( m=l ) * 

II-2 

11,403  (1) 

1.07 

Beaded  plate 

-b 

96 

19 

94  (1) 

1.02 

Flanged  plate 
(simple 
support ) 

-c 

862  (1) 

9 

907  (1) 

0.95 

Flanged  plate 
( clamped) 

-d 

1,840.(1) 

9 

1,853  (1) 

0.99 

Flanged  plate 

-e 

19,600  (1) 

II-2 

18,899  (1) 

1.04 

Lipped 

channel 

-f 

1,030 

22 

1,138  (10) 

0.91 

Lipped 

Z-section 

-g 

1,030 

22 

1,137  (10) 

0.91 

*  m  =  axial  half-wave  number 


3.  Geometry  of  Boeing  Test  Specimens  of  Composite  Plates 

In  Section  (7,  B)  of  the  main  text,  the  results  of  correlations  between 
tests  of  Boeing  specimens  of  composite  plates  and  the  analytical  predictions 
were  given  in  Table  7.3.  The  geometry  of  these  specimens  are  given  here 
in  Table  II. 2. 
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Table  II. 2  Composite  Plate  Thickness  Data  (Boeing) 


Plate 

Specimen 

Layer 
Thickness 
and  Total 
Thickness 

Layer 

Thicknesses 

( inch) 

Specimen 

(-D 

Specimen 

(-2) 

Specimen 

(-3) 

1  -  1 

Titanium 

.0195 

.0195 

.0195 

8A 

Adhesive 

.0185 

.0195 

.0125 

Composite 

.0650 

.0650 

.0650 

Total 

.10  3 

.104 

.097 

Titanium 

.0317 

.0317 

.0317 

8B 

Adhesive 

.0166 

.0176 

.0156 

Composite 

.0557 

.0557 

.0557 

Total 

.104 

.105 

.103 

Titanium 

.0340 

.0340 

.0340 

8C 

Adhesive 

.020 

.019 

.019 

Composite 

.0470 

.0470 

.0470 

Total 

.101 

.100 

.100 

Titanium 

.0504 

.0504 

.0504 

8D 

Adhesive 

.0171 

.0181 

.0181 

Composite 

.0375 

.0375 

.0375 

Total 

.105 

.106 

.106 

Composite 

.0332 

.0332 

.0332 

Adhesive 

.0107 

.0122 

.0122 

8K 

Titanium 

.0232 

.0232 

.0232 

Adhesive 

.0107 

.0122 

.0122 

Composite 

.0332 

.0332 

.0332 

Total 

.111 

.114 

.Ilk 

Composite 

.0274 

.0274 

.0274 

Adhesive 

.01825 

.01625 

.01575 

Of 

Titanium 

.0327 

.0327 

.0327 

Adhesive 

.01825 

.01625 

.01575 

Composite 

.0274 

.0274 

.0274 

Total 

.124 

.120 

.119 

Composite 

.0232 

.0232 

.0232 

Adhesive 

.02055 

.01955 

.01605 

8G 

Titanium 

.0395 

.0395 

.0395 

Adhesive 

.02055 

.01955 

.01605 

Composite 

.0232 

.0232 

.0232 

Total 

.127 

.125 

.118 
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Table  II-2  Continued 


— 

Composite 

.0168  .0168  .0168 

Adhesi ve 

.0127  .0147  .0132 

8H 

Titanium 

.0530  .0530  .0530 

Adhesive 

.0127  .011*7  .0132 

Composite 

.0168  .0168  .0168 

Total 

.112  .116  .113 

4.  Determination  of  Buckling  Load  by  Southwell  Plot 

As  mentioned  in  the  main  text,  in  some  cases  of  the  buckling  test  of 
composite  plate  specimens  which  results  are  given  in  Table  7.3,  it  is 
difficult  to  determine  the  buckling  load  from  a  load-deflection  curve. 

Some  of  the  buckling  loads  are  obtained  by  the  Southwell  plot.  Two  of 
such  plots  are  shown  in  Figures  II. 3  and  II. 4  for  reference.  In  Figure 
II. 3,  the  Southwell  plot  helped  to  decide  that  buckling  occurred  at  a 
higher  large -deflect ion  point  rather  than  at  a  lower  value  in  the  longitu¬ 
dinal  load-deflection  curve.  In  Figure  II. 4,  the  plot  helped  to  determine 
the  buckling  load  in  the  longitudinal  load-deflection  curve  which  seemed 
to  have  no  noticeable  buckling  point. 


References  for  the  Appendix 

II-l  Timoshenko,  S.  P. ,  and  Gere,  J.  M. ,  "Theory  of  Elastic  Stability," 
2nd  Ed.,  1961,  McGraw-Hill  Book  Company. 

II-2  Goodman,  S.,  "Elastic  Buckling  of  Outstanding  Flanges  Clamped  at  OKe 
Edge  and  Reinforced  by  Bulbs  at  the  Other,"  TN-1958,  October  1949, 
NACA,  Table  2,  p.  8. 
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Fig.  II. 1  Buckling  of  angle  section  (Ref.  II-l) 


rhirve  A:  Theoretical  curve  from  Ref.  II-l  for  buckling 

of  section  treated  as  two  separate  plate  elements 
and  each  one  is  simply  supported  at  three  sides 
and  the  fourth  side  free. 

Curve  B:  Theoretical  curve  from  Ref.  II-l  for  buckling 
of  section  treated  as  an  Euler  column  with  the 
ends  hinged. 

0  From  the  present  analysis 
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(-a)  7.50  0.400  2.80  0.200  Clamped 

'-b)  4.1+8  0.064  1.28  0.032  S.S. 


(-c)  3.20  Clamped 

(-d)  2.56  S.S. 


Beaded  plate 


Flanged  plate 


(-e)  Flanged  plate  a  =  10.0  (-f)  Lipped  channel  a  =  15.0 

(-g)  Lipped  Zee-section  a  =  15.0 


Fig.  II. 2  Geometry  of  sections  with  isotropic  beads  and  lips 

a  -  length  of  specimen 


•  Ml 


SOUTHWELL  PLOT 


Fig.  II -3  Longitudinal  and  lateral  load- do fleet ion  currea  of  flat  plate  teat  apeciwwn 
SO  No.  2  (Claap-claap- free- free  ease)  and  the  Southwell  plot 


Fig.  II-**  Longitudinal  and  lateral  load-deflection  currea  of  flat  plate  teat 

apeeiaen  SA  No.  1  (Claap-claiap-siinple-ainple  case)  and  the  Southwell 
plot 
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